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Abstract 

The muon anomalous magnetic moment and the hadronic vacuum polarization are 
examined using data analyzed within the framework of a suitably broken HLS model. The 
analysis relies on all available scan data samples and leaves provisionally aside the ex- 
isting ISR data. Our HLS model based global fit approach allows for a better check of 
consistency between data sets and we investigate how results depend on different strate- 
gies which may be followed. Relying on global fit qualities, we find several acceptable 
solutions leading to ambiguities in the reconstructed value for (a^th- Among these, 
the most conservative solution is ajf d,LO [HLS improved] = 687.72(4.63) x 10~ 10 and 
(o/i)th = 11659175.37(5.31) x 10 -10 corresponding to a 4.1<r significance for the dif- 
ference Aa /t = (a^) exp — (a^)th- It is also shown that the various contributions acces- 
sible through the model yield uniformly a factor 2 improvement of their uncertainty. The 
breaking procedure implemented in the HLS model is an extension of the former proce- 
dure based on a mechanism defined by Bando, Kugo and Yamawaki (BKY) . This yields 
a quite satisfactory simultaneous description of most e + e~ annihilation channels up to 
and including the (j> meson (7r + 7r~, 7r°7, 777, 7r + 7r~7r°, K + K~ , K°K°) and of a set of 
10 (mostly radiative) decay widths of light mesons. It also allows to achieve the proof 
of consistency between the e + e~ — > ■k + it~ annihilation and the r ± — > 7r ± 7r ^ decay 
and gives a solution to the reported problem concerning the measured partial width ratio 
T(0 -> K+K-)/T(<j) K°K ). Prospects for improving the VMD based estimates of 
a M are emphasized. 
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1 Introduction 



The muon anomalous magnetic moment is a physics piece of information which has 
been measured with the remarkable accuracy of 6.3 x 10 -10 [[UEl- From a theoretical point of 
view, a M is the sum of several contributions; the most prominent contributions can be predicted 
with a very high accuracy by the Standard Model. This covers the QED contribution which 
presently reaches an accuracy better than 1.6 x 10 -12 [3] or the electroweak contribution where 
the precision is now 1.8 x 10 -11 flU. The light-by-light contribution to a M is more complicated 
to estimate and is currently known with an accuracy of 2.6 x 10~ 10 ||5l . 

Another important contribution to a M is the hadronic vacuum polarization (HVP). Perturba- 
tive QCD allows to compute a part of this with an accuracy of the order 10 -11 ; this covers the 
high energy tail and the perturbative window between the J ftp and T resonance regions. For 
the region below this threshold, one is in the non-perturbative region of QCD where estimates 
of the hadronic VP cannot so far be directly derived from QCD, relying on first principles only. 
However, this may change in a future. Indeed, some recent progress in Lattice QCD [|6l|71[H 
gives hope that reliable calculations of the HVP are now in reach in the next years. They would 
be an important complement to the standard approaches, as well as to the approach presented 
here. 

One is, therefore, left with estimates numerically derived from experimental data. Indeed, 
it has been proved long ago that the contribution of an intermediate hadronic state Hi to a M is 
related with the annihilation cross section cth^s) = cr(e + e~ — > Hi) by : 



where K(s) is a known kernel [4] enhancing the weight of the low s region, close to the 
threshold of the final state Hi. Then, the total non-pertubative HVP can be estimated by 
a^{H) = a ti(Hi), where the sum extends over all final states Hi which can be reached in 
e + e~ annihilations. 

The accuracy of a^(Hi) is, of course, tightly related with the accuracy of the experimental 
data set used to perform numerically the integration shown above. When different data sets 
are available for a given annihilation channel Hi, a combination of the corresponding a M (ifj)'s 
is performed by weighting each estimate with the reported uncertainties affecting each data 
set, using standard statistical methods (see [J9!], for instance). Possible mismatches between the 
various estimates are accounted for by methods like the S-factor technics of the Particle Data 
Group ifTOl . In this approach, of course, the accuracy of each a^(Hi) is solely determined by 
all the measurements covering the channel Hi only, without any regard to the other channels 



This method succeeds in providing very precise values for the relevant contributions. Sum- 
ming up the non-perturbative HVP estimated this way with the rest, one obtains an estimate of 
quite comparable to the BNL average measurement [Q~l|2). However, the prediction based on 
e + e annihilation data or r decay data [fTTl IT2l H~3l IT4l IT31 l4l ITSTl exhibits a long-standing dis- 
crepancy; the exact value of this discrepancy has gone several times back and forth, depending 
on whether one trusts the r data based analyses or the scan e + e~ annihilation data, which are 
obviously more directly related with a^(H). With the advent of the high statistics data samples 




Hi (j + i). 
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collected using the Initial State Radiation (ISR) method lfT71[T8l[T9l , a precise value for this - 
possible - discrepancy has become harder to define unambiguously. 

In order to get a firm conclusion concerning the numerical difference between the measured 
and calculated values of the muon anomalous magnetic moment Aa^ = (a^) exp — (a^) th , one 
should first understand why r based and e + e~ based analyses differ; one should also understand 
the differences between scan data and ISR data and possibly the differences between the various 
available ISR data samples, as the KLOE samples fT7llT9ll and the BaBar sample [18] seem to 
lead to somewhat conflicting results. 

Anyway, while all proposed values for (a^th differ from the average for (a^) exp , the the- 
oretical uncertainties start to be comparable to the experimental one. Therefore, it becomes 
interesting to look for a method able to reduce the uncertainty on (a^) t h by simply using the 
existing data. It is also an important issue to have a framework where the properties of each 
data set can be examined. 

In order to cover the low energy regime of strong interactions, the most common approach 
is to use effective Lagrangians which preserve the symmetry properties of QCD. At very low 
energies, Chiral Perturbation Theory (ChPT) represents such a framework. However, the realm 
covered by the usual ChPT is very limited (not much greater than the r] mass); Resonance 
Chiral Perturbation Theory (R%PT) permits to go much deeper inside the resonance region; it 
thus defines a framework suited to study the non-perturbative hadronic VP (HVP). 

It was soon recognized [|20l that the coupling constants occuring at order p 4 in ChPT were 
saturated by low lying meson resonances of various kinds (vector, axial, scalar, pseudoscalar) 
as soon as they can contribute. This emphasized the role of the fundamental vector meson nonet 
and confirmed the relevance of the Vector Meson Dominance (VMD) concept in low energy 
physics. Soon after, [21] proved that the Hidden Local Symmetry (HLS) model 11221 1231 and 
the Resonance Chiral Perturbation Theory (R%PT) were equivalent. Therefore, one may think 
that the HLS model provides a convenient and constraining QCD inspired framework for an 
improved determination of the HVP. It is, therefore, quite legitimate to check wether the HLS 
model allows a better determination of the HVP than the usual method sketched above. 

The basic HLS model has an important limitation for HVP studies : The vector resonances 
entering the model are only those embodied in the lowest mass vector meson nonet. This cer- 
tainly limits upwards the relevant energy range to ~ 1.05 GeV, i.e. slightly above the 0(1020) 
meson mass; going beyond while staying within the standard HLS framework certainly entails 
uncontrolled uncertainties due to the contribution of higher mass vector meson nonets. 

However, relying on the standard method, one can estimate the contribution of the region 
y/s G [m n o, to 83.3% of the total HVP and show that its uncertainty is also a large fraction 
of the total HVP uncertainty : ~ 4 x 10~ 10 when using only scan data or ~ 2.7 x 10~ 10 when 
using also the recent ISR data samples. For comparison, the uncertainty provided by the region 
above ~ 1.05 GeV is ~ 4 x 10 -10 . Therefore, any significant improvement on the knowledge 
of {cL^th m the region \/s < 1.05 GeV is certainly valuable. 

The (basic) HLS model provides a framework where the interrelations between the various 
observed decay channels are made explicit. The point is that the use of an adequate model 
allows for a global fit strategy. All available cross-section data are used to constrain the model 
parameters, which in turn allows us to predict physical amplitudes. Therefore, if the model 
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provides a statistically acceptable common solution to some set H = {Hi} of different pro- 
cessesS, each covered by one or several data sets, the fit results can serve to reconstruct reliably 
the a p {Hi) (Hi eU). ' 

Indeed, if a global fit of the set % of the various data samples is successfully performed, then 
the parameter values and their error covariance matrix summarize reliably all the knowledge of 
the set "H, including the physics correlations among them. Then, all cross sections contained 
in H can be estimated with an information improved by having taken into account all the 
underlying physics correlations. 

With the present formulation, of the HLS model, the various a^(Hi) can be reliably and 
accurately determined up to ~ 1.05 GeV, just including the <p resonance region. All the rest 
should presently be estimated by the methods usual in this field. 

One can substantiate the benefits drawn for using such a global model : 

• As the model is global, it implies algebraic relations between the various channels it 
encompasses. Therefore, the accuracy of the estimate for a^(Hi) is determined by the 
statistics available for any channel Hi and also by the statistics associated with all the 
other channels contained in H,. 

For instance, the accuracy for a fJi (n + n^) is certainly determined by the available statis- 
tics for e + e~ — > but all other data, acting as constraints, also contribute to the 
accuracy for a M (?T + 7T~). This is the role of the e + e~ — > 777 or e + e~ — > tt°7 annihilation 
data, but also those of the decay width for — > 7/7 or of the dipion spectrum in the 
r — > ixixv decay, etc. 

Conversely, the accuracy for a M (?T 7), for instance, is not only governed by the statistics 
available for e + e~ — > 7r°7, but also by those provided by the e + e~ — > ix + n~ or e + e~ — > 
777 data, etc. 

Therefore, the improvement expected from a global model should affect simultaneously 
all the channels contained in H, and contributing to a M . 

• As the breaking procedure is global, it affects simultaneously all physics channels related 
with each other by the Lagrangian model. A successfull global fit thus implies that it is 
validated by the fit quality of the largest possible set of data samples. This high degree 
of consistency indicates that the breaking mode^l is not simply had hoc. 

• Any data set is certainly subject to specific systematics; however, taking into account 
that the study we plan relies on 45 different data sets covering 6 different annihilation 
channels, 10 partial width decays (taken from the Review of Particle Properties [10]) and 
some decay spectra^, one may consider the effects of correlated systematics reasonably 

'These can be cross sections or various kinds of meson partial widths, or also decay spectra. Indeed, any piece 
of information able to constrain the model parameters is valuable. 

2 We mean that the breaking procedure we define is certainly a model, but it is not intended to solve only one 
issue in isolation, like the consistency between e + e~ — > ir + ir~ and r — > ttttv, without any regard to the rest 
of the correlated physics. Stated otherwise, it is validated only if its consequences for the other related physics 
channels are accepted by the corresponding data. 

3 Actually, it affects the dipion spectrum in the decay r — > nirv and in the anomalous rj/rj' —> mry decays, 
among others. 
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well smeared out. Indeed, one may consider unlikely that the systematics affecting as 
many different objects can pile up. 

Basically, what is proposed is to introduce the theoretical prejudice represented by one 
formulation of the VMD concept in order to constrain the data beyond genuine statistical con- 
sistency of the various data samples referring to the same physics channel. It has already 
been shown [24] that theoretical (VMD) relationships among various channels are highly con- 
straining. The present work plans to better explore such a framework with a much improved 
modelling. 

Conceptually, the idea to include some theoretical prejudice in order to reduce the uncer- 
tainties on a M is not completely new. A method to complement the e + e~ — >■ n + n~ data with 
the constraints of analyticity, unitarity and chiral symmetry has been initiated by [|25l |26l |27l 
with the aim of improving the n + n~ contribution to a M , but this has not been finalized. 

For the present study, we have found appropriate to discard the data collected using the 
Initial State Radiation (ISR) method |fT7l[T8l[T9ll ; indeed, because of the complicated structure 
of their systematics, they almost certainly call for a more complicated statistical treatment than 
the usual e + e~ scan data. The use of ISR data will be addressed in a forthcoming publication. 

The HLS model ll22l l23l complemented with its anomalous sector 11281 provides a frame- 
work able to encompass a large realm of low energy physics. This anomalous sector will be 
referred to hereafter as FKTUY sector. The non-anomalous sector allows to cover most e + e~ 
annihilation channels and some r decays. Thanks to its anomalous sector, the same framework 
also includes the radiative decays of light flavor mesons with couplings of the forrrE VP'j and 
P77 and also several anomalous annihilation channels. Actually, up to the <p meson mass, the 
only identified channel which remains outside the scope of the HLS model is the e + e~ — > n tu 
annihilation channel, due to the large effect of high mass vector resonances [j29l [30) presently 
not included in the HLS model. 

However, in order to use the HLS model beyond rather qualitative studies and yield precise 
descriptions of experimental data, symmetry breaking procedures have to be implemented. A 
simple mechanism breaking the SU(3) flavor symmetry [|3D has been introduced, followed by 
several useful variants [1321 l33l l34l . Nonet symmetry breaking in the pseudoscalar sector has 
also been introduced by means of determinant terms 11331 . This breaking procedure has been 
shown to describe precisely the radiative decays of light mesons 11361 1371 and to meet [[381 all 
expectations of Chiral Perturbation Theory. 

In order to account for the reported mismatch between the pion form factor in e + e _ anni- 
hilation and in the r decay, it has been proposed [|3~9l to take into account loop effects. Indeed 
kaon loops produce a mixing of the neutral vector mesons which is a consequence for the 
K° — K ± mass splitting. These turn out to modify effectively the vector meson mass term by 
identified s-dependent terms. 

Introducing the physical vector fields which correspond to the eigenstates of the loop mod- 
ified vector meson mass matrix, provides a mixing mechanism of the triplet p° — cu — <fi system. 
In this change of fields the charged vector mesons remain unchanged. With this s-dependent 

4 In the following, we may denote by V and P any of respectively the vector or the pseudoscalar light flavor 
mesons. This does not rise ambiguities. 
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mixing of neutral vector mesons, the fit residuals to the pion form factor in e + e annihilations 
and in r decays did not exhibit any longer any mass dependence ll39~l ; thus this mechanism 



provides an important part of the solution to the so-called e + e~-r puzzlej. 

However, this solution is only partial. Indeed, if the dipion spectrum lineshape in the decay 
of the r lepton is clearly predicted [1391 1241 from e + e~ data, there is still some problem with 
its absolute magnitude. This issue has been found to be cured by allowing il a mass (5m) and 
a coupling (5g) difference between the neutral and charged p mesons, ii/ a rescaling of the r 
dipion spectra consistent with the reported uncertainties on the absolute scales of the various 
measured spectra [|40l HT1 @2). The results returned by fits did not lead to a significant mass 
difference^ but, instead, 8g and the fitted scales of the experimental spectra were found highly 
significant |[24|. 

However, the numerical values of these parameters (never more than a few percent) suggest 
that some unaccounted for isospin breaking effects have not yet been included. 

On the other hand, the HLS model supplemented with the SU(3)/U(3) breaking reminded 
above accounts successfully - and simultaneously - for the measured cross sections in the 
e + e~ — > 7r + 7r~, e + e~ — > 7r°7, e + e~ — > 777, e + e~ — > annihilation channels and for 

the additional set of 9 decay widths, especially the radiative decays of the form VP'j or -P77, 
needed in order to constrain more tightly the model. This has been proved in ll46l . However, 
as it stands, the HLS model fails to account for the annihilation channels e + e~ — > K + K~ 
and e + e~ — > K°K simultaneously. This is obviously related to the puzzling issue thoroughly 
discussed in ll47l concerning the branching fraction ratio <fi — > K + K~ /(f) — > K Tt°. The 
reported disagreement with theoretical expectations is found significant and amounts to a few 
percent. This also allows thinking that some isospin breaking effects are not yet fully accounted 



In the present paper, we define a symmetry breaking procedure which is nothing but an 
extension of the BKY mechanism referred to above, but including now breaking in the non- 
strange sectors. This mechanism is only an upgrade of the BKY mechanism and applies like- 
wise to the two different sectors (the so-called Ca and C v sectors) of the non-anomalous HLS 
Lagrangian. We show that the r scale issue is solved by breaking the Cy Lagrangian piece 
while the <fi — > K + K~ / <fi — > K°TC° puzzle yields its solution from applying the same mech- 
anism to the Ca Lagrangian piece. Stated otherwise, within the framework of the HLS model 
broken in this way, the e + e~-r and the </> — > KK puzzles appear as twin phenomena yielding 
parent explanations. 

Actually, equipped with this upgraded breaking mechanism, the HLS model provides a 
satisfactory description of all the physics information listed above, including now both e + e~ — > 
KK annihilations. 

Having discarded the 3 existing ISR data samples, a priori 45 different data sets of scan 
data are relevant for our present analysis. At each step of our analysis, we have checked the 
consistency of the various data samples with each other by relying, as strictly as possible, on 
the information provided by the various groups without any further assumption. We have found 

5 A similar result has been obtained in [16| relying rather on a p° — 7 mixing mechanism; it should be 
interesting to study a more general V — 7 mechanism supplementing the p° — cj — cf) mixing scheme. 

6 The mass difference following from fit corresponds to 8m — 0.68 ± 0.40 MeV is in accord with what is 
expected for the electromagnetic mass difference ll43]|441 of the p mesons 1451 dm ~ 0.81 MeV. 
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that 2 among them have a behavior not in agreement with what can be expected from the rest 
(43 data sets). One could have attempted to use them by weighting their contribution to the 
global x 2 ( a sort of S-factor); however, for now, we have preferred discarding them. Therefore 
our analysis relies on 43 data sets - mostly produced by the CMD-2 and SND Collaborations 
- and 10 accepted partial width information, which represents already an unusually large set of 
data consistently examined and satisfactorily understood. 

The paper is organized as follows. In Section [2] we briefly outine the basics of the HLS 
model and its various sectors. In Section[3] we define the upgraded breaking procedure which is 
a trivial extension to the u and d sectors of the BKY breaking scheme as redefined in If3"4l - the 
so-called "new scheme". Section [4] and Section \5\ examine the consequence for the modified 
BKY breaking scheme on the two different parts (Ca and Cy ) of the non-anomalous HLS 
Lagrangian. In Section [6] we first remind the loop mixing scheme [|39ll46l of the vector mesons 
and, next, construct the pion form factor in e + e" annihilation and in the decay of the r lepton. 
The condition F n (0) = 1 has some consequences for how parametrizing the Breit-Wigner 
amplitudes should be done for narrow objects like the to and mesons. Other topics are also 
examined : the direct conn coupling and the — > KK couplings. The anomalous sector is 
examined in SectionfT] where we also provide the expressions for the e + e~ — > 7r°7, e + e~ — » 777 
and e + e _ — > cross sections. The expression for the various couplings of the form 

P77 and VP'j are also derived; these are important ingredients for the set of radiative decays 
included into the HLS framework. 

We have found it appropriate to summarize the main features of the HLS model under the 
upgraded breaking scheme which underlies the present study; this is the matter of Section [8] 
Section [9] is devoted to list the different data sets available for each physics channel; in this 
Section, our fitting method, previously defined and used in [|39ll46ll24l . is reminded. 

At this point, we are in position to confront our model and the data. Section [TOlexamines the 
fit properties of the available e + e~ — > data and Section (TUreports on the simultaneous 

analysis of the e + e~ — > K + K~ and e + e~ — > K°K° annihilation data. The analysis of the 
e + e~ — > KK channels allows us to show how the problem raised by both — > KK decay 
widths is solved within the new release of the broken HLS model. 

Section d2l provides our analysis of the dipion spectrum in the r decay in conjunction with 
all e + e~ data. It is shown therein that e + e~ data and r data are fully reconciled; the precise 
mechanism solving this issue, somewhat unexpected, is exhibited. 

The short Section [13] is devoted to examining the exact structure of the lotxtx coupling and 
compare with similar results of other authors [ff8l [491. Similarly, another short Section [14] 
examines in some detail some properties of the 7r° — 77 — 77' mixing; it is shown here that 
the conclusions derived in ll38l about the mixing angles 9 and 9 8 introduced by ll50l I5T1 are 
confirmed, together with their relationship with the traditional singlet-octet mixing angle 9p. 
In Section [T5] one examines the fitted values of the parameters involved in the absolute scale 
of the FKTUY anomalous Lagrangian pieces and compare with existing estimates; this leads 
to the conclusion that the usual assumption C3 = C4 is consistent with data. 

Section |T6]is devoted to study in detail the consequences for our HLS model determination 
of the non-perturbative part of the photon hadronic vacuum polarization. This is found to yield 
much reduced uncertainties compared to estimates derived by the direct averaging of data. 

The consequence for g — 2 are also examined with the conclusion that the theoretical pre- 
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diction differs from the BNL measurement [2]. The significance of this difference is shown 
to stay in between 4.07a and 4.33o\ This looks an important improvement, as we are still not 
using the ISR data. 

Finally Section [FT] provides a summary of our conclusions and the perspectives. A large 
part of the formulae have been pushed inside several Appendices in order to ease as much as 
possible the reading of the main text. 



2 The HLS Lagrangian 

The Hidden Local Symmetry Model (HLS) has been presented in full detail in Il22ll and, 
more recently, in [23] . One can also find brief accounts in [|52~ll3"4l . 

Beside its non-anomalous sector, which allows to address most e + e~ annihilation chan- 
nels and some r decays up to about the <p meson mass 113911461 . the HLS Model also contains 
an anomalous (FKTUY) sector [28J which provides couplings of the form VVP, VPPP, 
^PPPyP^ or P77 among light flavor mesons. These are the key in order to incorporate 
within the HLS framework the radiative decays of the form VPj or P — > 77, or decays im- 
portantly influenced by the box anomaly like 77/77' — > 7r + 7r~7 (see 11531 [371 for instance). It has 
been shown that, while implementing (U(l)) nonet symmetry and SU(3) symmetry breakings, 
one reaches a remarkable agreement with data fl36l 1371 . 

The anomalous pieces of the HLS Model are also the key tool when dealing with annihi- 
lation processes like e + e~ — > 7r°7, e + e~ — > 777 or e + e~ — > as successfully shown in 

In order to be self-contained, and without going into unnecessary detail, let us briefly re- 
mind the salient features of the HLS Model relevant for the present purpose. 
One defines the £ fields by : 

£,r,l = exp [ia/ f a ] exp [MP J f n ] (1) 



where the scalar field a is usually eliminated by means of a suitable gauge choice [22J (the 
so-called unitary gauge). However the decay constant f a still survives through the ratio a = 
fl I fl which is a basic (free) ingredient of the HLS Model. The standard VMD Lagrangian 
corresponds to having a = 2. The pseudoscalar field matrix P : 



/ TT 3 + 7S" 8 + 7S m i " + 1 K+ ^ 

12, 1 

v K ~ r "V 3* + vi* y 

(2) 

contains singlet (P ) and octet (P 8 ) terms. By 7r 3 we denote the bare neutral pion field; the 
traditional naming 7r° will be devoted to the fully renormalized neutral pion field. On the other 
hand, the usual 77 and 7/ meson fields are (essentially) combinations of the r] 8 and r] fields 
shown in Eq. ©. 
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The HLS Lagrangian is defined by 



C-HLS 



La + a£ 



v 



-&Tt[(D,JUtl - D^ R ) 2 ] 



f 2 

~Tr[L — R] 2 



(3) 



t ^,21 



f 2 

T 



Tr[L + Rf 



where the covariant derivatives are given by : 

D^ L = d^ L - igV^ L + i£ L C^ 
D^r = d^ R - igV^ R + i^RKfj, 



(4) 



with 



= eQA, 



92 



cos9 w 

92 



-{T z - sin 2 6 W )Z, + ^{W;T + + W~T_) 



■ sin 2 OwZfj, 



(5) 



cos 9 W 

exhibiting the Z, W ± boson fields together with the photon field A^. The vector field matrix is 
given by : 

/ (p' + u 1 )/^ p + K* + \ 



V 



1 

V2 



V 



p 



K* 



-p 1 + u 1 



K* 
2 K*° 



(6) 



i 7 J 



The quark charge matrix Q is standard and the matrix T + = [T_]t is constructed out of 
matrix elements of the Cabibbo-Kobayashi-Maskawa matrix [|23l [39ll . One should note that 
the neutral charge entries of the vector field matrix V are expressed in terms of the so-called 
ideal fields (p 1 , uo 1 and (j) 1 ). 

In the expressions above, one observes the electric charge e, the universal vector coupling 



g and the weak coupling g 2 (related with the Fermi constant by g 2 = 2m w \j GfV^)- As the 
influence of the Z boson field is quite negligible in the physics we address, the Weinberg angle 
9 W plays no role. 

We do not present here the anomalous sectors which can be found in the original HLS 
literature [|22] |28] |23] . A summarized version, well suited to the present purpose, can be found 
in ll46l ) and will not be repeated. 

The non-anomalous Lagrangian Chls at lowest order in field derivatives can be found 
expanded in Il34ll52ll . Its r sector is explicitly given in [|39ll46l . 

The HLS Lagrangian fulfills a U(N f ) x U(N f ) symmetry rather than SU(N f ) x SU(N f ). 
The additional axial U(l) symmetry has several undesirable features Il35ll54l . especially a ninth 
light pseudoscalar meson. This symmetry can easily be reduced by adding appropriate terms to 
the effective Lagrangian. Defining ll22l the chiral field U = £|£r — ex P 2«P//^, this reduction 
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is achieved by adding determinant terms[35] to the HLS Lagrangian. After this operation, one 
gets [f3l : 

u 2 1 

£ = C-HLS + CtHooft = C-hls + y?7o + ^ Ad^o^o (7) 

where [i has obviously a mass dimension and A is dimensionless. In the following, the ad- 
ditional Lagrangian piece will not be modified while breaking symmetries. Actually, in the 
present work, one is only concerned by the perturbation of the pseudoscalar meson kinetic 
energy. 



3 The BKY-BOC Breaking of the HLS Lagrangian 

The HLS Lagrangian above is certainly an interesting and attractive framework. However, 
without introducing suitable mechanisms for symmetry breaking effects, one cannot account 
for the experimental data at the level of precision required by their accuracy. There is no 
unique way to implement such a mechanism within the HLS model and, actually, several SU(3) 
breaking schemes exist. The basic SU(3) symmetry breaking scheme has been proposed by 
Bando, Kugo and Yamawaki (BKY) 1 3 1 ] . It has, however, some undesirable properties which 
have motivated its modification. A first acceptable modification has been proposed by Bramon, 
Grau and Pancheri 11321 l33l and another one in [1341 . where the various schemes have been 
critically examined. Following this study, we prefer using in the following the so-called "new 
scheme" variant defined in [|3~4l : when referring to the BKY mechanism throughout this paper, 
we always mean the "new scheme" variant just mentioned. It will be referred to as either BKY 
or BKY-BOC. 

This breaking mechanism (BKY-BOC) has been examined in detail and its predictions - 
relying on fits to experimental data - have been found to meet the corresponding ChPT expec- 
tations ll38l at first order in the breaking parameters. It has also been extensively used in several 
successful studies performed on radiative decays of light mesons[36| and on e + e~ annihilation 
cross sections [|46l . Up to now, the BKY mechanism was limited to SU(3) symmetry breaking 
effects; the issue now is to examine its extension to isospin symmetry breaking. 

Briefly speaking, our variant of the BKY mechanism 041 turns out to define the broken 
non-anomalous HLS Lagrangian pieces by : 

C A = -^Tr[(L-R)X A ] 2 

f 2 (8) 
C v =-hL T r[(L + R)X v ] 2 

where Xa and Xy are matrices carrying the SU(3) symmetry breaking associated with, respec- 
tively, Ca and Cy. These are written as : 

X A = Diag(l, 1, z A ) 

(9) 

X v = Diag(l, l,Zy) 
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and departures of z A and zy from 1 account for SU(3) symmetry breaking effects in the La 
and Cy Lagrangian pieces^. A priori, these two parameters are unrelated and should be treated 
as independent of each other. 

In order to extend to isospin symmetry breaking, we propose to generalize Eqs. ® above 

to : 

( X A = Biag(q A ,y A ,z A ) 

(10) 

[ X v = Diag(q v ,y v ,z v ) . 

As isospin symmetry breaking is expected milder than SU(3) breaking, the additional 
breaking parameters are obviously expected to fulfill : 

\q A - 1\, \y A - 1| « \za - 1| , \q v - l\, \y v - 1\ « \z v - 1\ (11) 

In previous fits, performed with only SU(3) symmetry breaking, we got (see for instance [|39l 
l46l ) \z A — 1|, \zy — 1| ~ 0.5. Such ways of extending the BKY breaking mechanism have 
been already proposed within similar contexts [|55l |56l . 
We find appropriate to define : 

_ 1 , JU _ 1 , + A A /V 

q A /v - 1 + e A/v - H 

(12) 

_i , J _ i , S A/V - A A /y 
y A /v - 1 + e A/v - 1 H ^ 

exhibiting the sum and difference of e A j v and e A / v . Indeed, the expressions for most physical 
couplings are simpler in terms of these rather than in terms of e A / v and e A i v . 

As clear from Equations dHJ, the BKY breaking of the HLS Lagrangian exhibits a global 
character. It does not correspond to some systematic way of including specific breaking terms 
of given kind or order as done within ChPT. As the numerical values of the breaking parame- 
ters are phenomenologically derived from fits to a large set of experimental data, they account 
globally for several effects of different origin without any way to disentangle the various con- 
tributions. This remark is especially relevant for the breaking parameters corresponding to the 
u and d entries of X A and Xy which are small enough that several competing effects can mix 
togethe|| because of their relatively large magnitude, the SU(3) breaking effects can be more 
easily identified [[381 [36H . 

7 In the following X A and Xy are named breaking matrices; this convenient naming should not hide that the 
true breaking matrices are rather X A — 1 and Xy — 1. 

8 This may include effects due to the quark mass breaking and to electromagnetic corrections. It may also 
absorb corrections to hadronic vertices which can hardly be derived from an effective Lagrangian. 
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4 Breaking the Ca Lagrangian Piece 

The pseudoscalar kinetic energy term of the full (broken) Lagrangian is carried by C A + 
CtHooft- In terms of bare fields, it is : 

C K in = qAVAd-n + ■ <9vr~ + qA ^ yA dir 3 ■ dir 3 + q A z A dK + ■ OK' + y A z A dK° ■ dK° 



+ 



3 4 12 



6 2 

2 „,2 J2 „,2 



dr] ■ dr] (13) 



6 y/U V6 

which is clearly non-canonical. In order to restore the canonical structure, one should perform 
a change of fields. This is done in two steps, as in ll38l . 



4.1 First Step PS Field Renormalization 

The first step renormalization turns out to define the (step one) renormalized pseudoscalar 
field matrix P Rl in term of the bare one P by : 

P R i = X^PX 1 / 2 (14) 

The charged pion and both kaon terms in this expression actually undergo their final renormal- 
ization already at this (first) step. Indeed, the axial currents are defined as in ll38l and are given 
by: 

j; = -2U[Ti[T a X A d,PX A ] + \5 afi d^} (15) 

in terms of the bare fields and of the Gell-Mann matrices T a , normalized such as Tr[T a T b ] = 
5 a b /2. The ti ± , K ± and K°/TC° decay constants are defined by the relevant axial current 
matrix elements closed on the renormalized PS meson fields : < 0| ^ K \^% >= ifn/KQn- As 
one chooses the renormalized (charged) pion decay constant to coincide with its experimental 
central value [1571 (f w ± = 92.42 MeV), this turns out to impose q A y A = 1. At leading order in 
breaking parameters, this implies Y, A = 0. Then, the breaking matrix X A writes : 

X A = Diag(l + ^,l-^,z A ) (16) 

depending on only two free parameters (A^ and z A ). On the other hand, the kaon decay 
constant is : 

f K ± = v^I(i + (iv) 



One thus yields a marginal change compared to the previous BKY breaking scheme [1341 
(dealing only with SU(3) symmetry) as one got [/x±//vr±] 2 = z A . Anticipating somewhat on 
our numerical results, let us mention that the fits always return ~ (5. -r 6.)%, much larger 
than expected from solely an effect of the light quark mass difference ll58l : this will be further 
discussed in Subsection II 1.21 
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As clear from Eq. (fTol . the Xa matrix resembles the usual quark mass breaking matrix. 
However, the entry za is essentially related with the ratio [fx/ f-n] 2 — 1-5, while the corre- 
sponding entry in the quark mass breaking matrix is numerically ~ 20. Therefore, the corre- 
spondence between these two matrices is only formal. 

The following relationship defines some bare PS fields in terms of their (fully) renormalized 
partners : 



7T~ 



7T 



R 



K° 



'I 



4 



)K 



R 



(18) 



K 



R 



This produces changes going in opposite directions for the couplings involving the 
physical and K° mesons compared to their bare partners. This has a clear influence on 
the cross section ratio cr(e + e~ — > K + K~)/a(e + e~ — > K°K°). On the other hand, one also 
gets at leading order in breaking parameters the following relationship between some bare PS 
fields and the (first step) renormalized PS fields : 



7T 3 



71 



Ri 



A 



A 



A 



A „R 



7To 



y/2z A -l R , 



y/6 

A a B lz A + 2 R 

7=^% H Vr 

2^ 3 3 z A 18 



(19) 



V2z A -l R 



4.2 Second Step PS Field Renormalization 

While propagating the field redefinition displayed into Eqs. ([T8l[T9l in the expression for 
Cxin (Eq. (TT3l)). one should neglect second (and higher) order terms in the breaking parameters 
Aa and A. Indeed, both of these are expected small (of the order of a few percent at most); 
instead, as \z A — 1| is rather large (za — 1.5), we do not proceed alike with the SU(3) breaking 
term. Doing so, in terms of the (redefined) fields, the only surviving non-canonical piece £ ,8 
writes : 



2£ , 8 = [d^f + [drit? + ^ [(2 + y ) drff 1 + V2 (l 



1 

za 



(20) 



and is independent of A^. As we get - at leading order - the same dependence as before [|38l . 
the diagonalization procedure for this term is known (see Section 3.1 in lT38l0 . Let us only 
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recall the results in the present notation set 



where : 



Ri 



TT 



TT 



R 



1 + v cos 2 (3 R v sin /3 cos /3 R 



1 + u 



1 + v 



f sin <9 cos /3 R 1 + 1> sin 2 /3 R 

-^8 + — r-; w 



l + t> 



1 + V 



cos /3 



2^ + 1 
3(24 + 1) 



sin 



V2(z A - 1) 



3(24 + i; 



\ 



1 + A 



(24 + 1) 
34 



^A(24 + l) 



34 



(21) 



(22) 



It thus looks more appropriate to use v (~ A/2) rather than [|38l A as a breaking parameter, 
as it allows to work with simpler expressions, v is the first parameter in our model which 
exhibits the intricacy between U(l) and SU(3) breakings (A and za). The canonical PS fields - 
denoted by the superscript R - are finally defined by Eqs. (I18I19I21I) . 



4.3 The 7T° - 77 - rf Mixing 

The physically observed 77 and rf are traditionally described as mixtures of the singlet and 
octet PS fields i] 8 and 77 involving one mixing angle named here 9 p. Some authors, following 
Il59ll54ll prefer now using mixtures of the uu + dd and ss wave functions. However, as these 
two approaches are equivalent, we prefer sticking to the traditional description. 

Since ll5Tll50ll . it is admitted that the most appropriate ChPT description of the i] — 7]' mixing 
involves two decay constants (F° and F 8 ) and two mixing angles (0 O and 6*8). However, [|38l 
has shown how, within VMD, the usual octet-singlet mixing scheme connects with this new 
approach. In this reference, it was also shown that, relying on experimental data, the broken 
HLS model favors 9 = with a very good accuracy; this led to a relation between 9$ and 9p 
numerically close to 9 8 = 29 p. Comparing accepted ChPT numerical values for # 8 - like those 
in |j5TTl - with the one derived from 9p (determined in VMD fits) was found quite satisfactory. 
Moreover, it was shown that fits to experimental data lead to an algebraic relation of the form 
9 P = /(A, za). We will check whether this relation for 9 P still fits within the present form of 
our broken HLS model. 

As in all previous studies in this series, one could have limited oneself to considering only 
the ri — r]' mixing, decoupling this from the tt° sector. However, it is a classical matter of Chiral 
Perturbation Theory to address the issue of tt° — i] mixing, as this is related with the (light) 
quark mass difference [|60l . Therefore, it may look interesting to see if such a phenomenon 
could be exhibited from the experimental data we deal with. In this case, there is no reason not 
to address the issue of the relevance of a full tt° — 77 — rj mixing mechanism. We choose to 
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parametrize this PS mixing using [1611 : 

4 = - e 77 - e ' V 

< = cos#p(r7 + e7r ) + sin# P (r/ + e' tt°) (23) 

k Vr = - sin 6» P (?7 + e vr°) + cos P (»/ + e' tt°) 

where 7r^, 77^ and rf R are the already redefined fields (see Eqs. (|2TT i above), the physically 
observable mesons being 7r°, 77 and 77'. In the smooth limit of vanishing e and e', one recovers 
the usual 77 — 7/ mixing pattern with one (0p) mixing angle, while the pion field decouples. Even 
if one does not expect a large influence of e and e' in the full data set collection we consider, it 
does not harm to examine their effects and, if relevant impose e = e' = to the model. 

Finally, at leading order in breaking parameters, the pseudoscalar meson kinetic energy 
term is canonical when expressed in terms of the fully renormalized fields (those carrying a R 
subscript). 



4.4 About The 8 , 9 and 9 P Mixing Angles 

Let us define : 

' g*°( V ) = (l-2v) (-1-20%) 
v (2z A + V 2 



9» 
9 8 (v) 



1 - - 



3 (2z 



2^ v 



3 (24 + 1) l " 



(~ 1 - 10%) 
- 0.3%) 



(24) 



These functions tend to unity when the U A (1) symmetry is restored (A = 0). The property 
g 8 (v) ~ 1 is the simplest way to justify the approximation done in our previous works to 
parametrize nonet symmetry breaking by the parameter x (see, for instance, the discussion in 
ll38l ). Using these functions, one can derive from Eqs (fT3T > the following axial currents : 



J? = U \d^ R + A A g"(v) 



2^« + i=« 



Jf = U {^,4 + ZJ ^9» « - V2 Z -^g\v) (25) 



jv° 



0= d y R - v/2^fV(.) « + 2ZA + 1 - -* 



3 u v , ^ 3 
The mixing angles 9 8 , 9 [|5Tll50l yield the following expressions 



: 9 8 (v) d^l 



< 0\d>*J* \rf > 
tan 9 s = ; J| Q „ T « ,. : = tan (6 P - A) 



tan 8 



< 0\d»J*\r) > 

< 0|^JJ|t7 > 
~< 0\dfMrf > 



tan A = V2 



z A - 1 g°(v) 
2z A + lg8{v) 



tan 



B) 



, tani^v/2^M^ 



(26) 



z A + 2g\v) 
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One can easily check that g 8 (v)/g°(v) ~ 1 — A/2. A property to check from fits using the 
present form of the model is whether 9 is still consistent with zero ll38l . In this case, 9p is still 
no longer a free parameter, but fully determined by A and za, i-e. by breaking parameters and 
6 P tends to zero when the symmetries are restored. 

One should also note that the usual ChPT definition of the n° decay constant (< 0| J^ 3 |7r° >= 
iq^f^o) provides f^o = f w ±, not influenced by our isospin breaking procedure. However, as 
will be seen shortly - and as already emphasized in [|38l for the decays 77/7/ — >■ 77 - this is not 
the quantity actually involved in the decay n° — > 77. 



5 Breaking the Cy Lagrangian Piece 

The Cy Lagrangian, is defined by Eqs. (|8l-fl"2l. It yields the following vector meson mass 
term (m 2 = ag 2 f 2 ) : 



2 

_ m? 

I— mass 2 



1 + E v )p 2 + (1 + H v )uj + 2A vPl ■ u! + z V (fyj + 2(1 + £y)p + • p- (27) 



while keeping only the leading terms in the breaking parameters Ey and Ay (the K* mass term 
has been dropped out). As can be seen, the canonical structure of the mass term is broken by a 
p ■ uj term. 

In order to restore the canonical form of the mass term, one should perform a field re- 
definition in only the (p — u) sector. Interestingly, the requested transform is not a rotation 
but : 



Pi \ / PRi 



A 



v 



( h v u Rl 
\ (1 - h v )p Rl 



(28) 



makes the work when non-leading terms in Sy and Ay are neglected |^. In terms of the Ri 
renormalized fields, one gets : 

,2 



r - _ 

*-mass 2 



;i + Sy)p^ + (1 + Z v )u Rl + z v cj> Rl + 2(1 + Sy)p+ • p-\ (29) 
having renamed for convenience (pj = 4> Rl . A few remarks are worth being done : 



(1) Beside the two breaking parameters Sy and Ay, one gets a third free parameter hy 
which governs the mixing of p Rl and u Rl . 

Therefore, the exact content of isospin 1 ( Pl ) inside co Rl and of isospin (uj) inside p Rl 
should be extracted from data. 

(2) The masses for p Rl , u Rl and p 1 * 1 remain degenerated at leading order in the breaking 
parameters as the needed R\ change of fields results in a vector meson mass term inde- 
pendent of Ay. 

Therefore, even if one may legitimately think that breaking isospin symmetry inside Cy 
could result into a non zero (Lagrangian) mass difference [|24l Sm 2 between the charged 



9 Eq. d27l i can be diagonalized by a 45° rotation, however, this solution is physically unacceptable as it has not 
the requested smooth limit when Ay — > 0. 
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and the neutral p mesons, our breaking procedure rules out such a possibility at leading 
order in breaking parameters. Actually, electromagnetic corrections 03] [44], presently 
neglected, generate such a mass difference. Such a term has been considered in ||24~1 
but found numerically insignificant; preliminary studies within the present framework 
leading to the same conclusion, we have given up considering explicitly a p° — p ± mass 
difference. 

(3) The field transform (|29l ) propagates to the vector meson kinetic energy by generating 
a term of the form Aydp^du)^ which breaks the canonical structure of the kinetic 
energy. This is a classical issue [62 1 known to imply the occurence of wave-function 
renormalization factors [|62l which are absorbed into the effective couplings defined by 
the Lagrangian vertices. In our case, they are certainly absorbed in our breaking parame- 
ters. This is exactly the same issue which arises in the electroweak Standard Model with 
the 7 — Z° mixing. This has been investigated in detail within Z° lineshape studies (at 
the one plus two-loop level) and by the LEP experiments. The same issue also appears 
when treating the 7 — p° mixing and has been discussed in [16]. 

The second step renormalization of the vector meson fields, which accounts for loop effects 
lT3~9lf46l, is considered below. 



6 The Fully Broken Non-Anomalous HLS Lagrangian 



For definiteness, we name (abusively) from now on "HLS Lagrangian" the full expression 
given in Eq. CD), i.e. including the determinant terms. The HLS Lagrangian is explicitly 
provided in Appendix[A] dropping out for conciseness all terms not relevant for the purpose of 
the present study. 

At the present step - which does not still include the (second step) redefinition of the neutral 
vector fields [|39l [461 - several remarks are worth being done : 

• The vector meson masses occuring in the Lagrangian fulfill m 2 = m 2 ± = m^. Thus, no 
mass splitting is generated, except for the meson. 

• The couplings pnir undergo isospin breaking (Ey) but remain strictly identical for the 
charged and neutral p mesons. Instead, a direct cunn coupling is generated; it is propor- 
tional to (1 — h v )A v . 

• The p° — 7 and p^ — W ± transition amplitudes^ may slightly differ, as /iyAy/3 should 
not exceed a few percent level. 

Therefore, non-vanishing 5m 2 = m 2 — m 2 ± and 5g p7Tn = g p o n +^- — g p ±- K ±- K o, as stated in 
E4l . are not derived by extending the X^/Xy breaking scheme to include isospin symmetry 



Therefore, non-vanishing 5m 2 and Sg pn7r are not the way followed by the (broken) HLS 
model in order to account for the (slightly) different normalizations of the pion form factor in 
t decays and in e + e~ annihilations. The actual mechanism at work is emphasized below. 

'"Compare / P7 and f p \y as given by Eqs. d93l ) and d94l i. 

11 As stated above, electromagnetic corrections contribute to generate a non-vanishing 5m 2 without, however, 
a significant influence on the fit properties. 




16 



6.1 Loop Mixing of Vector Mesons 



As remarked in [39J, pseudoscalar loops modify the vector mass matrix by s-dependent 
terms. In this way, the p, to and <p squared masses become s-dependent through contributions 
at real s of analytic functions^ namely the KK loops and, for the p, also the charged pion 
loop. Conceptually, this turns out to remark that the inverse vector meson propagator written 



n(s) in order to exhibit the loop effects, can be thought as D v (s 



s — mr 



[s), reflecting the running character of the vector meson squared mass. 



More important, however, is that this s-dependent mass matrix becomes non-diagonal, 
showing that, at one-loop order, the p, u and <f> (corresponding here to the Ri renormalized vec- 
tor fields) are no-longer mass eigenstates. Mass eigenstates can easily be constructed by stan- 
dard perturbative methods [|63l as shown in [|39l ; one observes that they become s-dependent. 

This mass matrix can be written : 



M 2 {s) 



M 2 i 



s) + 5M 2 (s) 



where0 : 



M 2 (s) 



Diagfm 2 + II, 



(30) 



(3D 



is treated as the unperturbed part of the squared mass matrix. The pion loop is weighted by 
the square of the p^nn coupling constant (see Eq.(l92l) in Appendix A) and has been included 
in the p Rl entry as Il 7r7r (s) is not really small in the timelike region. Instead, as the uj Rl im 
coupling is first order in Ay, the pion loop contribution to the co^ entry should be neglected 
(~ 0(A\)). The values for these (Higgs-Kibble) masses can be found in Eq. (1931) ; they fulfill 



On the other hand, 5M 2 (s) is given by : 



5M 2 (s) 



^-p4> i £u4> i £<j> 



(32) 



and contains only the perturbations generated by kaon loop effects. The kaon loop transition 
from a vector meson V to another one V has been denoted eyv ■ 

One should note |[39l that M 2 is an analytic function of s satisfying the (so-called) hermi- 
tian analyticity condition : M 2 (s*) = [M 2 (s)] t . 

The entries of these matrices are appropriately parametrized in terms of : 



n 



n (s) 



e2 (s)= n+(s) + n (s) 



g 



pTTTT 



9pn 



(33) 



12 Actually, the anomalous FKTUY Lagrangian and the Yang-Mills terms contribute respectively with VP and 
VV loops; one can consider their influence absorbed in the subtraction polynomials of the PP loops l39ll . 
13 Entries are ordered respectively pn lt oj^ and ^>r 1 . 
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where U.'(s) denotes the amputated pion loop, while 11+ (s) and ITq(s) denote, respectively, 
the amputated charged and neutral kaon loops; their analytic expressions can be found in the 
Appendices of 11391 . ei(s) ^(s) do not contain symmetry breaking terms beyond the effects of 
the kaon mass splitting. The expressions for the entries in SM 2 (s) are given in Appendix B 
and show this dependence explicitly (see Eqs. (197T)). 

One can construct, as in ||3~9~l , the eigensolutions of M 2 . These are the final (step two) 
renormalized vector fields denoted respectively by p R , uo R and <p R and are related with their Ri 
partners by : 

/ PRi \ ( PR- auj R + fi<pR \ 



V 



+ ap R + 

\<Pr- PpR - l^R J 
where the s-dependent mixing angles are defined by : 

e 



(34) 



a{s) 
7(s) 



-pUJ 



Xp — X u 



e p<$> 



Xp — Xq 



(35) 



Xu> — Xqs 

using the eigenvalues of M 2 (at first order) : 

X p (s) = m 2 + TL^(s) + e p (s) , A w (s) = ml + ej^s) , X^s) = + e^(s) (36) 

The e p (s), e w (s) and e^(s) quantities, defined in Eqs.(l97l), only depend on the kaon loop func- 
tions and on breaking parameters. 



6.2 The Pion Form Factor in e + e Annihilations and in r Decays 

The pion form factor in the ^ decay to tx ± t\^v t can easily be derived from the Lagrangian 
piece C T given in Eq. (l94l) : 



where : 



FKs) 



f( 1 + s ^z^ 



(37) 



f fas) = /; - n w ( s ) 



D p (s) = s-ml-UL(s 



2/1 iv "\ ™2 _ „ J f2/ 



(38) 
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and the loop functions are 



IMs) 



r ag, 



(1 - |(1 + E v ))4(s) + , -d + S r ))/. /v (.s) 



24 



Pw(s) 



24 



P P (s) 



(39) 

where ^(s) and £k(s) denote respectively the amputated charged pion and kaon loops, P w (s) 
and P p (s) being subtraction polynomials. In order to fulfill current conservation, these poly- 
nomials should vanish at s = 0. Here, as in former studies [|3~9ll46ll24l . identifying P ± P T and 
P ± P° loops has been found numerically justified. 

If one compares with the corresponding formulae in [24] (Subsection 2.1.1), one sees that 
Sm 2 and Sg - supposed to reflect different properties of the charged and neutral p mesons - 
have been deleted. As the loop functions vanish at s = 0, one clearly has F£(0) = 1. 

The pion form factor in e + e~ annihilations is not as simply derived. One needs first to 
propagate the transformation in Eq. (1341 into the Lagrangian Eq. (|92l and collect all contribu- 
tions to, respectively, p R , co R and <p R . In this way, the V — 7 couplings associated with the fully 
renormalized vector fields become : 



-z v , 



P 



f 



3 

agfl. 



1 + E v + 3(1 - h v )A v - 3a(s) + V2^{s)z v ) 
-V2z v + 3/3(s) + 7 (s)) 



(40) 



including the mixing angle contributions. Using the Lagrangian pieces given in Eqs.(|99l), one 
can construct easily the pion form factor : 



where : 



2 K v 3 



(41) 



gpim — -^k J- + 



f [(1-MAy 



ag at \ 



The loop corrected V — 7 transitions amplitudes F v (s) are defined by : 



(42) 



= fv ~ ny 7 (s) , (V = p% u R , <p R ) 



(43) 



with the s-dependent loop terms IL/ 7 (s) being defined in Appendix C. All ily 7 (s) are re- 
quested to vanish at s = because of current conservation. The inverse p propagator D p (s) is 
defined by (see Eq. (1361) ) : 



D p o(s) 



s - X p (s) 



m. 



(44) 
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F*{s) (1=1) 


Non-Resonant Term 


[l-|(l + Ev)] 




p Mass Squared 


ag 2 fl{l + E v ) 


a«? 2 / 2 (l + E y ) 


7T7T Coupling 




f(l + 5V) 


Amplitudes /; & f; 


a^ 2 (l + E v ) 




J P 

f T 
J P 


1+ 3 + 3 + 3 " y 


Renormalization factor of 
KK couplings 


1 

za 





Table 1 : Comparison of the pion form factor information in r decay and in e + e~ annihilation. 
Second column lists only isospin 1 related information. In the last entry of the rightmost data 
column, the upper sign refers to K + K~ pairs, the lower to K°K°. 
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As the p self-mass H pp (s) vanishes at s = 0, one certainly has D p o(0) = —m p . Con- 
cerning the to and mesons, one can correspondingly write their inverse propagators as : 

D v {s) = s-m 2 v - U vv (s) , (ml = m 2 p , mj = z v m 2 p ) (45) 

and one can legitimately assume their self-energies to also vanish at s = 0. Then, D^O) = 
—m 2 and A/>(0) = — zym 2 should certainly be fulfilled. However (most of) the u self-energy 
cannot be computed in closed form and the 3-pion part of <p self-energy too. Therefore, con- 
venient forms for their propagators should be considered. This issue is readdressed just below 
for both mesons. 

At this step, it is of concern to compare the properties of the isospin 1 part of F*(s) with 
F£(s) . The most important pieces of information are listed in TableQ] The difference displayed 
for the non-resonant term is tiny. One can see that there is no mass difference between the 
charged and neutral p mesons, nor different couplings to a pion pair. Instead, most of the 
difference is actually carried out by the transition amplitudes (see the fifth data line in Tabled)) 
which are significantly s-dependent, as can be inferred from Figures 6 and 7 in [46 j . 

Finally, it is interesting to note that the renormalization factor introduced in couplings in- 
volving a kaon pair plays in opposite directions for charged and neutral kaon pairs. 

6.3 The ujtttt Direct Coupling and the Condition F^(0) = 1 

As can be seen from Eqs. (|42]> . the fully broken HLS model reveals a total coupling of the 
to to a pion pair given by : 

9uwx = y K 1 ~ h V )A V - a(s)} . 

This expression illustrates that the cunn coupling in our model is a priori a superposition of 
a direct isospin breaking term and of another piece generated by vector meson mixing through 
kaon loops. This kind of sharing has been emphasized several times [48, 49J. The full data set 
we use should give the most precise and motivated estimate for these two pieces as this is still 
presently controversial [|48ll49ll64l . 

On the other hand, the parametrization of the cu contribution to the pion form factor may 
pose a conceptual problem related with the condition F*(0) = 1 which is worth addressing. 

The pseudoscalar meson loops which enter the W transition amplitudes (see Eqs. (l32l) . 
(1351) and (1971) ) behave as 0(s) near the origin. The running vector meson masses (see Eqs. 
(l36l) ) are such that \ p (s) — A u ,(s) vanishes at the origin, while the two other differences which 
come into Eqs. (l36l) tend to a non-zero constant. Therefore, ab initio, the mixing angles are 
expected to fulfill : 

0(0) = 7 (0) = , «(0) ^ (46) 

Even if clear in the previous publications using the loop mixing mechanism (Figure 7 in 
||39l or Figure 6 in [[461 clearly show that a(0) ~ —5% ), this was not explicitly pointed out. 
Therefore, the s-dependent ujhh coupling generated by loop mixing^ does not vanish at the 
origin. This has some consequences. 

14 Actually, the non-identically vanishing e\{s) function providing the vector meson mixing via loops is also 
generated by isospin symmetry breaking, however in the pseudoscalar sector. 
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Indeed, using Eqs. (l40l) and (l42l) . and the vanishing properties of the functions IIy 7 (s), 
Uvv(s) and /3(s) at the origin, one gets : 



^ e (0) = 1+ g 

when keeping only the first-order terms in breaking parameters. 

As already discussed at the end of the previous Subsection, it is motivated to assume the u 
self-energy 11^ (s) vanishing at the origin. Moreover, this allows to stay consistent with the 
so-called "Node theorem" Il65ll66l . Then, the inverse propagator D u (s) = s — m 2 — n ww (s) 
fulfills -0^(0) = — m 2 . This provides the vanishing of the last bracket in the formula above 
and, thus, F*(0) = 1, whatever the values for hy, Ay and a(0). 

However, in most applications, for objects carrying such a narrow width as the u and <\> 
mesons, one generally uses approximate inverse propagators, e.g. eithei0 : 

D v (s) = s -m\ + im v f v (BW a ) or D v (s) = s - m v + iy/sty (BW b ). 

with values for my and Ty either taken from the Review of Particle Properties or extracted 
from one's fits. Then, with either of these Breit-Wigner lineshapes, the condition F%(0) = 1 
is not necessarily fulfilled. From our model results, this condition is even violated at a few 
percent level. However, it is easy to check that either of : 

D U] (s) = s-m 2 -(ml - m 2 - im u T w ) (BW a ') 

m p 1 

(remember that m 2 p = m^) and : 

D u (s) = s-m 2 p - -^{ml -m 2 - i^sY^) (BW b ') 
1 m z p p 

certainly cures this disease. This turns out to parametrize the u self-energy n ww (s) with an 
ansatz which satisfies its vanishing at the origin. 

It is worth stressing that using standard Breit-Wigner lineshapes or their modified partners 
provides practically unchanged fit results. This is due to the fact that the p and u masses (with 
tilde or not) are close to each other, and then, the factor s/m 2 p is very well approximated by 1 
all along the sensitive region of the to peak. 

In order to substantiate the possible changes, we have run our code using BW a and BW a ' as 
inverse to propagators. As a typical example of modification, one can compare m w = 782.44 ± 
0.06 MeV and f u = 8.46 ± 0.09 MeV while fitting with BW a , and m u = 782.49 ± 0.06 MeV 
and r w = 8.36 ± 0.08 MeV when using instead BW a '. For definiteness, in the fits presented in 
this paper, D u (s) will be modified as just explained. As /3(0) = 0, the pion form factor value 
at s = is not sensistive to how the <p propagator is approximated. 

15 Within the ongoing discussion, phenomenological values -e.g. not derived from the broken HLS model 
parameters values- for vector mesons masses and widths, are denoted with a tilde symbol in order to avoid 
confusion. 



a(0) 



h v )A v 



m r 



Au(O) 



(47) 
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Even if our choice is motivated, others are certainly possible as exemplified in 11131 14911 . 
Transposed to our model, the just mentioned choice would turn out to weight the full to con- 
tribution to the pion form factor by a factor s/m^ or s/m 2 p which restores F%(0) = 1. The 
behavior of this choice is identical to ours, basically because m w and m p are very close to each 
other. 



6.4 The Charged and Neutral Kaon Form Factors 

We give here the annihilation cross sections/form factors within the extended BKY-BOC 
breaking of the HLS Lagrangian. Cross sections and form factors are related through : 



a(e + e- ->• PP) 



— , ona 



m Q 3 p\F e P (s) 



3s 5 /2 



(48) 



for any meson pair PP. q P = ys — Amp /2 is the P momentum in the center-of-mass system. 
The kaon form factors are given by : 



FUs) 



-^-(2 + Zy + 2Sy + 2Ay - ^(2 + Zy) 

6za 2 



-^-(1 - Zy + £ y - Ay + -^(1 - Zy)) 

bZA 2 



9 p k+k-F p1 (s) 9ujk+k-F w1 (s) g<j,K+K~F^(s) 



D P (s) 



D u (s) D^s) 
3 pK ^F P1 {s) g^-oF^is) g^j^F^is) 



D P {s) 



D^s) 
(49) 

where the 7 — V transition amplitudes Fy 1 have been already defined (see Eq. (l43l). The 
VKK couplings can be read off from the corresponding Lagrangian pieces (Eqs. (1 1001) and 

dm). 

The kaon form factors fulfill : 



However, it is easy to check that these conditions are both fulfilled, only if : 



z v m 2 p 



A,(0) D^O) 



-1 



(50) 



(51) 



Therefore a fixed width Breit-Wigner shape for the should be adapted as already dis- 
cussed for the u. 



6.5 Parametrization of the 4> Propagator 

As for the pion form factor, in order to fulfill F^ (0) = 1 and Ft(0) = 0, one should 
impose that the u and inverse propagators at s = are equal in magnitude and opposite in 
sign to their respective Lagrangian masses (m| = zym^ = zym p ). Here again, this turns out 
to parametrize the full self-energy li^(s) by an ansatz vanishing at s = 0. For the two-body 
loops, this is well known 11391 ; however, the three-body loop is not known in closed form (as 
for IL^s)). 
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However, in contrast with the case for tu, using 



DJs) = s - z v m 2 -(ml - z v m 2 - im^T^) (BW a ') (52) 

for the 4> inverse propagator, instead of the usual (fixed width) BW. form, should be further 
commented, as ^fz~y~m p significantly differs from normally fitted^! (e.g. with BW a ). 

Even if anticipating on our fit results, it is worth discussing this matter right now. As 
far as cross sections are concerned, the two kinds of fits provide almost identical results. In 
order to yield this result, almost all parameters vary within errors except foi0 zy, which could 
have been expected. However, it will be shown that this change has a marginal influence on 
all information of physics importance. Anyway, such kind of information is interesting as it 
provides a hint on the model dependence of numerical results. Therefore, it has been of concern 
to compare results obtained with either of BW a and BW a ', when appropriate. 

Before closing this Section, one may note that, at the <\> peak location (y/s ~ 1020 Mev), 
the modified Breit-Wigner lineshape provides : 

z v m 2 + -^(m? - z v m 2 ) ~ [(1.020) MeV] 2 
1 ZymA 

which explains why the fit remains successful when using B W a ' and also why zy should change 
correspondingly, taking into account that m 2 p cannot much vary. The fit quality of the e + e~ — >■ 
KK cross sections will illustrate the validity of this parametrization of the <p propagator. 



6.6 The Coulomb Interaction Factor 

Beyond modelling, there is an important issue to discuss when dealing with the charged 
kaon form factor. In the decay <p — > K + K~, and more generally as close to the KK threshold, 
one has to take into account the Coulomb interaction among the emerging charged kaons. This 
has been first addressed in ll67l and recently readdressed (and corrected) in [PTTl . The net result 
of this effect is to multiply the charged kaon cross section by the Coulomb factoi0 : 



Z(s) 



1 + na e 



l + v 2 1 Js-4m 2 



2v 



K± 



(53) 



In [68J, and later in ll69l , the cross section for charged kaons is multiplied by Z(s)/Z(m^ ) ). 
This turns out to consider the Coulomb interaction as a breaking mechanism which affects the 
charged kaon sector and not the neutral one; as the corresponding <f> branching fractions are fit 
independently, this should not affect their results. One may just have to remark that this turns 
out to incorporate the Coulomb effects inside the corresponding estimates for the <p — > K + K~ 
branching fraction. 



16 More substantially, with appropriate fits, one yields ^fzyrrip ~ 925 MeV, while a direct fit yields ~ 1020 
MeV ! 

17 In fits with BW a for the meson, one gets zy = 1.368 ± 0.005, while with BW a ' the fit returns zy = 
1.472 ±0.001. 

18 Actually, the full electromagnetic correction factor is more complicated, but the main effect comes from 
the Coulomb factor. One assumes that the kaon data which have been submitted to fit have been appropriately 
corrected for soft photon corrections, which allows to cancel out the term named C, in [47 1. 
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6.7 About the -> K + R-/(J) -> Ratio 

Up to well defined phase space factors generated by the kaon mass splitting, the partial 
width ratio <fi — )■ K + K~ / cf) — >■ K 7C° is the square of the corresponding s-dependent effective 
coupling ratio. Neglecting for each coupling corrections terms of order greater than 1, one can 
derive from Eqs. (11001 ) and (11011 ) : 

= - i* , ~~'i~ 7 !'l 1 1 - A *i - - 1 1 - < 54 > 

where the last equation follows from remarking (see Figure 7 in [14610 that the mixing angle 
/3(s) - defined by Eq. (|35l) - is negligibly small compared to y2zy in the <p mass region. 
Therefore, this mechanism proposes a way for this ratio to depart from unity. 

In their throughout study of the — > K + K~ /<$ — > K°TC° ratio, the authors of [|4~71 exam- 
ined this issue using several other mechanisms than this one and concluded that none of them 
was able to accomodate a coupling constant ratio smaller than one (in absolute magnitude). 
The global fit, based on the suitably broken HLS model, provides a new approach. In this 
framework, the determination of is constrained by both e + e~ — > KK annihilation cross 
sections separately, and by some more light meson anomalous decays, which also depend on 



7 The HLS Anomalous Sector 

In order to treat radiative decays, i.e. the VP'j couplings, and some important annihila- 
tion channels (namely e + e~ — > 7r°7, e + e~ — > 777 and e + e~ — > within the HLS 
framework, one needs to incorporate the appropriate Lagrangian pieces. These are given by 
the Wess-Zumino-Witten (WZW) terms iTTOl 1711 which traditionally account for the triangle 
(AAP) and box (APPP) anomalies, together with the FKTUY Lagrangian pieces J28ll23l : 

4 

£-anom. = £-W ZW + 51 C *^« (^5) 
i=l 

where the four q are constants left unconstrained by theory [28]. A closer examination of the 
FKTUY Lagrangian allows to identify five different pieces - listed in Appendix O- and one 
then remarks that the accessible physics is sensitive to the difference c\ — c 2 and not to each of 
them separately. One is then left a priori with three unconstrained parameters [|23l . 

When no breaking is at work, the amplitudes for the couplings^ Po77 an d Po 7r+7r 7 at the 
chiral point - computed within the FKTUY-HLS frameworlo- coincide with those directly 
derived from the WZW piece in isolation [|46l . Due to a sign erroi0 in the FKTUY Lagrangian 
piece Cavp, it was asserted in [46] that the constraint c 3 = c 4 was mandatory in order to 

19 Here and in the following Pq denotes either of the tt°, r\ and r/' mesons. 

20 e.g. using Eq. (l55l l and the V — 7 transitions provided by the non-anomalous HLS Lagrangian. 
21 We gratefully ackowledge B. Kubis (HISKP, Bonn University) for having kindly pointed out the issue. 
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recover this property. Actually, this property is automatically satisfied Il22ll23l . In addition, we 
have verified that this property is maintained within our fully broken HLS model. 

However, the condition C3 = C4, which is fulfilled by VMD models [23j is successful and 
only turns out to reduce the freedom in fits. Nevertheless, one has examined relaxing this 
condition and found that our fit results are well compatible with the constraint c 3 = c 4 . 

7.1 Breaking the Anomalous HLS Lagrangian 

At this step, the anomalous HLS Lagrangian can be written : 



with pieces listed in Appendix |Dj As for the non-anomalous HLS Lagrangian, each among 
these pieces may undergo specific symmetry breaking independently of each other. This may 
lead to plenty of free parameters as illustrated by M. Hashimoto [56J who implemented com- 
bined SU(3) and Isospin symmetry breakings in the anomalous sector. 

A simpler mechanism has also been proposed for SU(3) breaking by Bramon, Grau and 
Pancheri ||3~2l [331 ; however, this was insufficient to account for both K*^' ^ — > K^' '^ de- 
cay widths. In 11361 l39l it was proposed to supplement it with a breaking of the vector field 
matrix resembling a vector field redefinition. Quite unexpectedly, this provides a (successful) 
parametrization for the K* radiative partial widths identical to those proposed by G. Morpurgo 
IT721 within a completely different context. Interestingly, this combined mechanism leaves to- 
tally unaffected the other sectors of the Cvvp piece we deal with; this is well accepted by all 
data considered ll36l l39l . This combined breaking mechanism has been studied in detail ll46l 
for all pieces of C anom , with similar conclusions. 

The combined breaking mechanism, as presented in ll46l . has been examined by combining 
SU(3) and Isospin symmetry breakings using the complete data set discussed below within the 
minimization code underlying the present study. It was concluded that possible Isospin symme- 
try breaking effects - not propagated from the field redefinitions provided by non-anomalous 
HLS Lagrangian breaking - provide invisible effects. It was then decided to neglect this addi- 
tional possible source of Isospin symmetry breaking, as the parameter freedom it gives is found 
useless. 

Therefore, for sake of clarity, one only quotes the specific forms for the decay amplitudes 
*[±,o] _i, Kl±,o]y^ referring the interested reader to [[461 for more information. 

As a summary, our dealing with the anomalous sector - except for the limited K* sector - 
involves only 3 parameters : c\ — c 2 and c 3 and c 4 ; former studies ll46l l24l |3~71 remain valid, as 
the condition C4 — C3 = is well accepted by the data, as will be shown shortly. 

7.2 Radiative Couplings 

For what concerns the radiative decays of light mesons and the e + e~ — > P7 annihilation 
processes, one needs Caap and an effective piece named C' AVP defined below. 

In terms of the final renormalized pseudoscalar fields and assuming the ir° — i] — rj mixing 
defined in Section |4l one can write : 



anom. 



£-AAP + £-APPP + £-VVP + C-VPPP 



(56) 



cm 



7T 



(1 - Cije^^daApd^Ay g^— + g my 



2V3 



n 




(57) 
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At leading order in breaking parameters, the coefficients gp ory are given b>@ : 



9n0 



77 



5Vy7 
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3 
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5Aa _e 
6 + v/3l 
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sin Op > + 
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13za(1- 
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l + 2z A 
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y/2- 



sm 
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(58) 



These clearly depend on the breaking parameters A a, z a and v (the PS nonet symmetry break- 
ing) and on the 7r° — 77 — 77' mixing scheme (see Eqs. (|24l) ). especially on the singlet-octet 
mixing angle Op. One should note that / 7r /(7 7r o 77 is another way to define the neutral pion decay 
constant. The other equations also illustrate that the so-called octet and singlet decay constants 
as derivable from there have little to do with the standardly defined ones, i.e. from the currents 
in Eqs. (|25l ). This question has raised some confusion which motivated the study in ||38l . 

In order to treat the e + e~ — > annihilation process the part of the C A ppp La- 

grangian describing the so-called box anomalies is needed. This can be written : 



E 
with : 



APPP 



J 7T 



4 (Cl 



- c 2 + c 4 



d a ir dpTT 



(59) 
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T]7T~*~7T 7 



S'?/7r+7r - 7 




1 + 2vz A 



V3 

1- ZA 




1 - VZ A 



AVP 



Eqs. (|58l) and (1601) show how the triangle and box anomaly amplitudes behave under 
isospin, SU(3) and PS nonet symmetry breakings. One should especially note the intricacy 
of SU(3) and PS nonet symmetry breakings. 

In order to derive the radiative decay couplings, an effective Lagrangian has been built up 
from Cyyp and the non-anomalous Lagrangian in the same way as in [46J. This can be written 
in terms of the renormalized Ri fields : 

(;<>""" i) n ,i) n .\., with g= 



GF, 



/IV 



P=ir°, r], rf 



gppj d^ 1 + g Pu)J d^ 1 + g P<h d, 



Ri 



(61) 



22 One could expand the (1 + v) 1 factor and keep only the contributions of orders 1 and v. However, in the 
present case, it does not simplify the expressions. 
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The expression for the various coupling constants gyp 1 can be found in Appendix El In order 
to derive the physical couplings, one should first apply the transformation given in Eq.(l34T) and 
then collect the various contributions to each of the (neutral) pr, lor and 4>r. 

Concerning the AVP couplings, it is quite interesting to compare the expressions in Eqs. 
(1 1 081 ) — (1 1 1 Ob with the corresponding ones in f|39ll38"l , derived using an approximate expression 
for nonet symmetry breaking^ (the x parameter in the quoted papers). Indeed, the three vari- 
ants by which nonet breaking occurs (see Eqs. (|107l l) are close together and can reasonably well 
approximated by x e ff = 1 — v ~ I — A/2. 

7.3 Breaking the VVP and VPPP Anomalous Lagrangians 

The VPPP anomalous Lagrangian is given by : 

' C VPPP = -%D<r°» { [g° pn d u 7r° + g^rj + +g° pv ,d u r ] '} pf 



+ 

with D 



3g(ci - C2 - c 3 ) 

J 7T 



(62) 



where one has limited oneself to display the V P n + n sector. The leading terms of the cou- 
plings occuring in this expression are given in Appendix |F] 



The Cvvp Lagrangian piece plays an important role in the annihilation process e + e — > 
7T°7r + 7r~. Its relevant part is : 
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(63) 
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2^fl 



(64) 



~A 



When going from ii^-renormalized to the fully renormalized vector fields R, one has to take 
some care with attributing the s-dependence between the two neutral fields of each monomial 
in the second line of Eq. (|63l) . This should be tracked for each R\ field while applying Eq. 
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In order to restore the condition C3 7^ C4, one should simply make in [39] the replacement C3 — > (03 + Ci)/2. 
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7.4 The e + e — > Po7 Annihilation Cross Sections 

Using the Lagrangian pieces given above, the transition amplitudes 7* — > P7 can be written 
similarly to 11461 : 



T( 7 * -> P o7 ) = iY 



9 ° 3 t° A Kp { s ) ~ (1 - c 4)£p 



e MI/a/ V,(gK^(p) , P = tt , 77, 7/ 

(65) 

where Y = —a em N c / irf^ has been factored out. q is the incoming photon momentum (q 2 = s), 
p the outgoing photon momentum (p 2 = 0) and N c = 3. The pieces provided by Caap are0 : 



9tt° 
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(66) 



using the gp o77 couplings defined in Eqs. (|58~T) . where the (1 — c 4 ) has been factored out. The 
resonance contributions are gathered in K Po (s) : 



K Po (s) 
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Vi=p R ,LO R ,4> R 



H^(s)F v ^(s) 
DvAs) 



P = tt , m 7/ 



(67) 



where the Hy° (s) - given in AppendixlGl- are the resonance couplings to P 7 an d the F v r Js) 
are the V — 7 transition amplitudes defined in Eq. (|43T ). The Dy^s) are the vector mesons 
inverse propagators already encountered. The cross sections can then be written : 



a(e + e- P o7 ) - 



Po 



I^Po» 



(68) 



7.5 The e + e — > 7r°7r + 7r Annihilation Cross Section 

Following as closely as before the notations in [|46l . the amplitude for the 7* — > is 
given by : 



T ( 7 * ^ Tr+rf) = [T sym ( s ) + T ferfc (s) + T AVP (s)] ^e^p^pp 



(69) 



where £^(9) (g 2 = s) is the (heavy) photon polarization vector. T sym is the symmetric part 
of the amplitude (in terms of the pn 'final' states), while Tb r k (denoted T p in [46]) breaks this 
symmetry. We have found appropriate to introduce separately the contribution Tavp(s) to 
the full amplitude generated by the Cavp Lagrangian piece (see Eqs. (11051 )); its first term is 



24 The corresponding expressions given in [46 1 carry a missprint : Each of the right-hand sides of Eqs. (41) is 
missing a factor of 2. 
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symmetric in terms of the pit 'final' states. One ha 



25 . 
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where all parameters and functions have been already defined, except for the iVj(s) functions 
which are given and commented in Appendix|H] One has kept as much as possible the notations 
used in ll46l in order to exhibit the effects of our additional isospin symmetry breaking effects 
by simple inspection. Finally, Tavp(s) identically vanishes when c 4 = c 3 . 
The differential cross section writes : 



d 2 cr(e + e — > 7T + 7r 7T°) a 



em 2 



dx dy 1927T 



s G(x } y)\T sym (s) + T hrk (s) + T AVP (s) 



(71) 



using the (x and y) parametrization proposed by E. Kuraev and Z. Siligadze [1731 who provided 
the kernel function G(x, y) reminded in Appendix IHl Note also that each of T sym (s), T brk (s) 
and Tavp(s) also depend on x and y. 



8 Ugraded Breaking of the HLS Model : A Summary 

In the former studies performed along the present lines [|39] @6l|24]|, roughly speaking, one 
incorporated nonet symmetry and SU(3) symmetry breaking in the pseudoscalar (PS) sector. 
In the vector meson sector, only SU(3) symmetry breaking was considered. 

However, some important effects can be already attributed to isospin breaking effects in the 
PS sector. Indeed, it is the non-vanishing character of the mixing "angles" a(s) and f3(s) which 
induces s-dependent p — cu and p — <\> mixings at the one loop level. This non-vanishing of the 
a(s) and /3(s) functions proceeds from the kaon mass splitting which breaks the symmetry be- 
tween the neutral and charged kaon loops and, then, allows to choose the analytic function e\ (s) 
as non-identically vanishing. Therefore, except for the u — <p system which would mix anyway 
at one loop, the full loop mixing mechanism for vector mesons is the prominent consequence 
for this limited account of isospin breaking^. 

This quite limited breaking scheme, allows already for a good account [|39l l46l l24l of the 
available data. However, within the realm accessible to the HLS model, two experimental 
issues remain unsolved : 

25 The contribution was wrongly omitted in the study ll46l ; the error was due to having missed that the two 

occurences of the function 7 in the numerator in the last Equation (II 15b come with two different arguments (s_| 

and s). The authors of the study [46 1 apologize for this inconvenience. 

26 Actually, as noted in previous works [ 39"], VP and VV loops contribute to the vector meson mixing; the 
effect of these additional loops can be considered as absorbed by the subtraction polynomials of the kaon loops. 
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i) The dipion spectrum lineshape in r decays is consistent with expectations from e + e 
annihilations 113911461 . but not its absolute scale [|24l . 

ii) The partial width ratio T(0 — > K + K~)/T(<f) — > K°K°) is found inconsistent with all re- 
ported expectations [|4"7ll . Obviously, this inconsistency propagates to the corresponding 
e + e~ annihilation cross sections. 

The first topic has been shown to get a satisfactory - but not perfect - solution by allowing 
some difference between p° and p^ meson properties to be fitted from data. If the effect of 
a non-vanishing 5m? = m 2 p0 — m 2 p± was found small, those generated by a non-vanishing 
& 9 — 9 P °n+Tr- — 9p±n°TT± was found especially significant [24]. Moreover, some rescaling of the 
t spectra, consistent with the reported experimental scale uncertainties remained unavoidable. 

The second topic is experimentally addressed by considering ll68l l69l that the Coulomb 
interaction^ plays as a symmetry breaking mechanism which modifies the SU(3) relationship 

94>k+k- = g^Koj? between coupling constants to g^+K- = 3^ f» \l z { m l)- Thi s ap- 
proach, which turns out to consider the Coulomb interaction as some breaking effect, may look 
unsatisfactory; anyway, it does not fit with our breaking scheme. 

These two issues motivated an upgrade of the breaking scheme of the HLS model in order 
to check whether an acceptable solution can be derived. The extension to isospin breaking of 
the BKY-BOC breaking mechanism is a priori an obvious candidate to examine. This has been 
done in the preceding Sections with several interesting conclusions, which can be summarized 
as follows : 

j) One does not find any signal for a mass or a coupling difference between the p° and p^ 1 
mesons^! However the coupling difference between p — 7 and p — W might be enforced 
with respect to [139114611241 if the breaking parameter product /iy Ay is found significantly 
non-zero (see Table [D, 

jj) Everything goes as if the universal coupling g remains unchanged in the anomalous sec- 
tor, while one observes that g is effectively modified to g(l + Ey) for the whole non- 
anomalous sector. Therefore, isospin breaking in the HLS model generates some mild 
disconnection between anomalous and non-anomalous processes which needs to be ex- 
plored. 

jjj) The partial width ratio T(<f) — > K + K~)/T(<p — > K°K°) is found subject to isospin 
breaking in a novel way compared with the various possibilities examined in H71 . 

Topics j and jj are both important for scale issues. Indeed, by disconnecting somewhat 
more than before the ratio of transition amplitudes p — 7 and p — W, one allows the HLS model 
to get more freedom for the purpose to account for scale issues. More important, both r and 
e + e~ physics share the same universal coupling (g(l + £y)), but it is no longer common with 
the scale of the anomalous processes which remains governed by g. Moreover, none among the 
anomalous couplings, all displayed in several of the Appendices, exhibits a dependence upon 

27 The function Z(s) in Eq. d53l 

28 Electromagnetic effects beyond the HLS model and the BKY breaking scheme may, of course, change a little 
bit this picture; however, the phenomenological consequences of letting free this mass difference are known to be 
negligible [24] as reminded before. 



31 



Ey. Stated otherwise, the anomalous couplings - which fix the scales of the anomalous meson 
decay and annihilation processes - no longer constrain the non-anomalous process scales as 
sharply as formerly assumed [|39ll46ll24l . 

Concerning the topics ii and jjj, it should be stressed that the parameter governing the 
change of this ratio is not involved only in the ratio. Indeed, each of the e + e~ — > K + K~ and 
e + e~ — > K°K° cross sections should keep valid absolute scales separately. Moreover, as clear 
from Appendices |El |F1 iGl and IH1 and from Eqs. (|58T ) and (1601 given above, this change of scale 
should also fit with all anomalous processes, including the n° — > 77 partial width, now within 
the partial width data sample submitted to the global fit. 

Before ending up this Section and this Part, let us remark that the upgraded breaking of 
the HLS model allows to address the question of the n° — 77 — rf mixing in an unusually large 
context. Moreover, as seen in Subsection l6.3[ the exact structure of the cunn coupling discussed 
several times in the literature [@8] |64] |49) can also be examined within the largest possible data 
set. 

A last remark is worth being emphasized. The scale treatment and the partial width ratio 
quoted in i and ii, within the upgraded breaking of the HLS model show up as two different 
aspects of the same mechanism. Indeed, the former proceeds from applying the extended 
BKY-BOC breaking scheme to C v , while the latter follows from applying the same mechanism 
to C A . 

9 The Data Sets and Their Handling 

In this Section, we outline the data sets submitted to the global fit and the way correlated and 
uncorrelated uncertainties are dealt with. Nothing really new is involved here compared to what 
is already stated in ||39l |46l 1241 . except for the data sets associated with the e + e~ — > K + K~ 
and e + e~ — > K°K° cross sections. One may, thus, consider that this Section is, to a large 
extent, a simple reminder provided in order to ease the reading of the present paper. 

9.1 The e e -> tt + tt + Data 

Four data sets have been collected recently in Novosibirsk at the VEPP2M ring. The first 
one [ 7411751 . covering the region from about 600 to 960 MeV, is claimed to carry a remarquably 
small systematic error (0.6%). Later, CMD-2 has published two additional data sets, one 
11761 - covering the energy region from 600 to 970 MeV - is supposed to reach a systematic 
error of 0.8%, and a second set [1771 closer to the threshold region (from 370 to 520 MeV) 
has an estimated systematic error of 0.7%. On the other hand, the SND collaboration has 
published ItTBl a data set covering the invariant mass region from 370 to 970 MeV. Except for 
the two data points closest to threshold which carry a sizable systematic error (3.2%), a reported 
systematic uncertainty of 1.3% affects this spectrum. These four data sets may be referred to 
in the following as "new timelike data" 11391 . 

When dealing with these data sets, statistical and uncorrelated systematic uncertainties have 
been added in quadrature as usual. However, these four data sets also carry a common corre- 
lated systematic uncertainty estimated to 0.4% which affects all of them in the same way [1791 . 
This is accounted for by modifying appropriately the covariance matrix as outlined in ||39l l46l 
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- see also Subsection 19.71 below - and by accounting for the data set to data set correlations. 
This is performed by treating these four data sets altogether, as if they were subsets of a single 
(merged) data set. 

In order to be complete, we have also included in our fit all data on the pion form factor 
collected formerly by the OLYA and CMD Collaborations as tabulated in (801 and the DM1 
data |f8Tfl collected at ACO (Orsay). These data will be referred to globally as "old timelike 
data". The systematic uncertainties carried by OLYA data (4%) and CMD (2%) contain an un- 
corrected part which has been added in quadrature to the reported statistical errors. A common 
correlated part of the systematics, conservatively estimated [|79l to 1%, has been dealt with ap- 
propriately. Instead, the accuracy of the DM1 data set being poor and its weight marginal, we 
did not find any need to go beyond the published uncorrected errors. 

9.2 The e + e- ->> (tt%)7 Data 

Since 1999, several data sets on the anomalous annihilation channels e + e~ — > 7r°7 and 
e + e~ — > 777 have been made available by the CMD-2 and SND Collaborations. In our analysis, 
we only use the provided data points up to y/s = 1.05 GeV. 

The first one used is the data set from CMD-2 11821 on the 777 final state (77 — > 7r + 7r~7r ) 
which carries a systematic error of 4.8%. CMD-2 has also provided [83J a second data set on 
the 777 final state, tagged with the decay mode 77 — > 3n°. The systematic uncertainty carried 
by this sample is estimated to 6.1% and 4.1% for, respectively, the energy regions below and 
above 950 MeV. More recently, CMD-2 has also published two more data sets [84] covering 
both the (7r°/^)7 final states, tagged with the 2-photon decay modes, in the energy region from 
600 to 1380 MeV. These are reported to carry a 6 % systematic error. 

The SND Collaboration has recently published ll8~5l two different data sets for the 777 final 
state with an estimated systematic uncertainty of ~ 4.8 %. The first one covers the energy 
region from 600 to 1360 MeV and the second from 755 to 1055 MeV. A sample covering the 
energy range from 600 to 970 MeV for the ir° — > 77 decay mode was also published ||86ll . 
Finally, two data sets for both (7c°/r])'~f final states with 14 data points (from 985 MeV to 1039 
MeV) from SND [87J are also available; these exhibit the much lower systematic error of 2.5%. 

Altogether, these two Collaborations have provided 86 measurement points for the e + e~ — > 
7r°7 cross section and 182 for e + e~ — > 777 for y^i < 1.05 GeV. Preliminary analyses [|46l did 
not reveal any need to split up correlated and uncorrected parts of the systematic errors for 
the (rj/n )^ data samples. Nevertheless, we have made a few checks by comparing fit results 
derived by adding in quadrature statistical and systematic uncertainties with fit results derived 
assuming the reported systematic error to be 100% bin-to-bin correlated. We did not observe 
any significant difference. Therefore, when analyzing the e + e~ — » (ir /^)^ data, the reported 
statistical and systematic uncertainties have been simply added in quadrature as in [|4~6l. 

9.3 The e+e- ->> 7r°7r + 7r Data 

This channel is important as it provides a single place where the box anomaly sector [17011711 
is present. Other physics channels involving the box anomaly in the 77/77' sectors exist (77/77' — > 
7r + 7r~7) and may be relevant. However, the overall experimental situation is unclear [13711461 . 
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even if the Crystal Barrel data sample [1531 may look secure. Therefore, we find preferable to 
wait for confirmation with new data samples which could come from BES and KLOE. 

There are several published data sets for the e + e~ — > annihilation channel with 

various statistical and systematic uncertainties. We first included in our data sample the data 
sets collected by CMD-2 which consist of a measured sample covering the to region [1751 
affected with a global scale uncertainty of 1.3% and two others which cover the region with 
a reported scale error of, respectively, 4.6% fl8~8l and 1.9% f|89l . The most recent CMD-2 data 
sample [|90ll also covers the <fi region with a scale uncertainty of 2.5%. 

SND has published two spectra covering altogether the region from 0.44 to 1.38 GeV, the 
former below 980 MeV fl9T1 . the latter above [92] . For both data samples, the correlated part of 
the systematic uncertainty has been extracted in order to be treated as a scale uncertainty (3.4 
% for [|9T| and 5% for [1921 , respectively); the uncorrelated parts have been added in quadrature 
with the reported statistical errors. 

Former data sets are also considered which cover the region in between the to and <p peaks 
where physics constraints are valuable. The most useful has been collected by the ND Collab- 
oration with 10% systematics and can be found in 11931 , the latter is a small data sample from 
CMD ||94l providing 5 measurement points with 15% systematics in the intermediate region. 
Concerning these two complementary data samples, we perform as in [|46ll and do not extract 
the correlated part of the systematics as the accuracy is poor enough that this could not lead 
to visible effects in global fits. Finally, there also exists a small data sample from DM1 [1931 
which has been used for illustrative purposes only [|46l. 

The analysis of these data samples has been performed in [1461 ; however, as the N 5 term 
which contributes to the cross section (see Eq. (f70"T )) was missing, the analysis is redone and 
the conclusions revisited. 



9.4 The r± -> t&tPv t Data 

In the collection of data samples submitted to global fitting, we also use the ALEPH ||40l , 
CLEO [|42l and BELLE [41] data sets. When dealing with r data, it is important to note that 
the relevant quantity, sensitive to the spectrum lineshape and to its absolute normalization is 
given by : 

1 dV^js) _ p 1 dN{s) 

T T ds ™N ds K ] 

where T T is the full r width, B n7r the branching ratio to ixnu, and l/NdN(s) / ds is the normal- 
ized spectrum of yields as measured by the various experiments. 

The data published by the ALEPH Collaboration correspond directly to the quantity shown 
in the left-hand side of Eq. d72l) . Instead, each of CLEO and BELLE has published separately 
the normalized spectrum of yields and the measured branching ratio B wn . In the r data han- 
dling, we have considered the reported uncertainties on these measured B^'s as bin-to-bin 
correlated scale uncertainties; these come into the various \ 2 associated with each data set in 
the way reminded in Subsection 19.71 Stated otherwise, they are no longer fitted as previously 
done B21. 

Following closely the experimental information provided by HOI . iHTTl . 11421 . the scale un- 
certainties have been estimated to 0.51% (ALEPH), 1.53% (Belle) and 1.74% (CLEO). On the 
other hand, a possible absolute energy scale uncertainty of 0.9% r.m.s. affecting the CLEO data 
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sample ll42ll has not been found significant [|39l l24l and is not considered in the present study. 
All these experiments have provided their statistical and systematic error covariance matrices; 
these are the main ingredient of the x 2 functions used in the fits. 

As the HLS model relies on the lowest mass vector meson nonet only, it cannot access 
T T which is therefore taken from the Review of Particle Properties [96|. Finally, our model 
provides [1391 : 

dT^Js) \V ud \ 



ds 
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and F£(s) is given in Eq. (1371) . Isospin symmetry breaking specific of the r decay will be 
considered and taken into account as emphasized in Section [T2l 

Of course, the published r spectra extend much beyond the validity range of the HLS 
model, as this presently stands. Therefore, when using it, we have to truncate at some s value. 
Consistency with the treatment of scan data would imply a truncation at 1 .05 GeV. However, 
various studies BT1 l24l showing the behavior of fit residual clearly observe that ALEPH data 
on the one hand and Belle and CLEO data, on the other hand, exhibit inconsistent behavior 
starting in the 0.9 4- 1. GeV region. Therefore, we have preferred truncating the spectrum at 1. 
GeV, where the three spectra are in reasonable agreement with each other. 



9.5 The e+e ->> KK Data 

Several data sets have been collected by the CMD-2 and SND Collaborations on both anni- 
hilation cross sections e + e" — > K + K~ and e + e~ — > K°K° . Here also, we have discarded the 
data points above 1.05 GeV. 

The oldest data sets, published by CMD-2 [88], provide the spectra for both the neutral 
and charged decay channels with a systematic uncertainty of 4%. Recently CMD-2 has rean- 
alyzed four data sets for the neutral decay mode [97] getting small systematic errors (1.7%). 
More recently, CMD-2 has also published two scans of the charged mode spectrum [|9~8l with a 
systematic uncertainty of 2.2%. 

On the other hand, SND has published in 2001 several data sets [[681 : 2 for the charged 
decay channel with a systematic error of 7. 1 %, 2 data sets in the neutral mode with K s — > 
and 2 more with K$ — > with respectively 4.2% and 4.0% systematics. 

The quoted systematics are treated as correlated scale uncertainty as outlined in Subsection 
l977lbelow. 



9.6 The Partial Width Data Set 

In order to work out the fit procedure and get enough constraints on the physics parameters 
of the model, an important input is the set of decay partial widths If39l . All decay modes of 
the form VP^f and P77 not related with the cross sections listed above should be considered. 



35 



This covers the radiative partial widths — > n ±, y, 7/ — > wy and — > 7/7 on the one hand 
and (7/ / 77/77°) — >■ 77 on the other hand. They have been extracted from the Review of Particle 
Properties [|96l . The accepted values for radiative partial widths for K*^ — > K ± j and K*° — >■ 
i^°7 have also to be used 11961 . 

As the currently available data on e + e~ — > n + ix~ stop slightly below 1 GeV, the phase 
of the (f) — > n + ix~ amplitude and its branching ratio as measured by SND [99] are relevant 
pieces of information, not included in the above listed annihilation dataEl In contrast, the 
corresponding information for the cu meson is irrelevant as it is fully contained in the amplitude 
for e + e~ — > n + ix~ (see Eq. (O) and is already part of the data sample. 

With respect to former studies within the same framework, the only new piece of informa- 
tion included in the fit data set is the partial width 7r° — > 77. Indeed, as can be seen from Eqs. 
(1581) . the corresponding amplitude may constrain as well as the e + e~ — > KK annihilation 
amplitudes. 

In fits involving all the above quoted annihilation channels, one has no longer to consider 
the leptonic widths (p°/w/0) — > e + e~ and the decay widths (p° /uj/4>) — > (rj/n )^ as they are 
essentially extracted from some of the cross sections listed above which permanently enter our 
fit procedure. 

Therefore, the additional decay information to be used as input to final fits represents in 
total 10 more pieces of information. 

9.7 Outline of the Fit Procedure (The Method) 

For all data sets listed above, one always has at one's disposal the statistical error covariance 
matrix. For scan data, this may include the uncorrected part of the systematic errors; if not 
done at start, enough information is generally provided to allow one to perform this (quadratic) 
sum. In the case of r data, the systematic error covariance matrix may be provided by the 
experimental groups (as ALEPH [40J, for instance). 

In this case, for each group of data sets (7r + 7r~, 7r°7, 777, 7r + 7r~7r°, K + K~, K°TC° , n ± TT u) 
one computes the partial x 2 '■ 

x \ = (771 - M) T V~\m - M) (Experiment # 1) (75) 

using matrix notations, and denoting by m and M the measurement vector and the correspond- 
ing model function vector. V is the error covariance matrix already referred to. The function 
to minimize is simply the sum of the xf- 

Actually, this is the procedure to estimate x] when the corresponding data sample is not 
subject to an overall scale uncertainty. If such a scale uncertainty takes place for some data set, 
one should perform a modification. 

Let us assume that the data set i is subject to a scale uncertainty; this is supposed^ to be a 
random variable e(0, a) of zero mean (unbiased) and with r.m.s. cr, independent of s. Then any 

29 However, one might have to be cautious with these data. Indeed, as emphasized in ll46l - see Section 13 
therein - the single piece of information truely model independent is the product BeeB^. Therefore separate 
values for B ee and B^ir, given as "experimental" values in the various releases of the Review of Particle Properties, 
are actually model dependent to an unknown extent. 

30 In practical use, a data set # i, subject to a scale uncertainty Aj,o is supposed to have been corrected in order 
to absorb a possible bias; this is the reason why the corresponding random variable is supposed unbiased, e.g. 
carrying zero mean. If not, Eq. (ITol should be modified by performing — > A; — A^o- 



36 



fit corresponds to getting one sampling of e(0, a), named Aj. In this case, Eq. (1751) should be 
modified to : 

X l = [m-M- AXif V- 1 [m-M- AX*] + -± (76) 

where [100] A is traditionally the vector of the model values M and the other notations are 
obvious. One can solve for A, which turns out to perform the change : 



V- 1 =^ W~Ha 



2\ 



V + a AA 



T 



'2 



v ~ l ~ i + ^v-W V ~ lA){V ~ lA)T <?7) 



in Eq. (|75T) . The modified covariance matrix W depends on the vector A. As just stated, the 
best motivated choice for the vector A is the model function A = M. However, this implies 
a recursive determination of the modified covariance matrix, and, therefore, recalculating (or 
inverting) large matrices at each step of the minimization procedure (several hundreds of times 
for each fit attempt). It happens, however, probably because the experimental data we deal 
with are already accurate enough, that choosing A = m (i.e. the measurement vector of the 
corresponding experiment) does not sensitively affect the results and strongly improves the 
convergence speed of the minimization procedure @6). Therefore, unless otherwise stated, we 
always perform this approximation. 



9.8 The Discarded Data Sets 

There exists data sets which have been discarded for the present study. The most important 
are the three data sets collected using the Initial State Radiation (ISR) method by the KLOE 
lfl7l [T9l and BaBar [18J Collaborations. These suppose a specific statistical treatment as the 
structure of the reported systematic errors is much more complex than for any set of scan data. 
The method used in ll46l for KLOE 2008 data ifTTll allows to deal with, but should be studied 
carefully with each ISR data set separately. 

In order to keep clear the message of the present study, we prefer avoiding using now data 
sets invoking delicate statistical methods. Therefore, the ISR data sets ifTTl [T9l [T8l will be 
treated in a forthcoming publication. Because of their high statistics, if well understood, these 
data samples may improve the physics results derived by using the model and the fit procedure 
presented in this study. 

Other data sets could have been useful : 

• Those providing the pion form factor in the spacelike region close to s = H101U1021 . 
Indeed such data could severely constrain the pion form factor in the threshold region. 
This was illustrated in [39] where an archaic form of our model has been used. However, 
we gave up using them - especially H1011I - because there is some suspicion concerning 
their estimated overall scale. Such a kind of data would nevertheless help in getting more 
precise information on g — 2. 

• More data involving the box anomaly, especially in the 77/77' sectors may also help in 
constraining the model parameters. For instance, the dipion spectra in 77/77' — > 7r + 7r~7 
provide such information. Some available data collected in 11371 , especially those for 
7/ — > 7r + 7r~7 provided by the Crystal Barrel Collaboration [53], might be considered 
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sometime. However, new data sets on this subject, with larger statistics and better sys- 
tematics should come from the KLOE and BES Collaborations, especially concerning 
the decay process 77 — > n + n~'~f. These are certainly more easy to handle than the 
e + e~ — > T]n + n^ annihilation data which in fine carry the same physics information. 

9.9 The Physics Parameter Set 

It looks appropriate to give the list of the free model parameters to be fitted from data. The 
model parameters are of various kinds : 

• The basic HLS (4) parameters : the universal vector coupling g; the relative weight a of 
the Lagrangian pieces Ca and £y, expected a ~ 2 from most VMD models; finally the 
weights c 3 , C4 and c\ — C2 of the anomalous FKTUY Lagrangian pieces to be added to the 
HLS Lagrangian in order to address the full set of data outlined in the above Subsections. 

• SU(3) breaking parameters which modifies the physics content of the HLS Lagrangian 
(za, zy and zt), together with the parameter named A which accounts for nonet symmetry 
breaking in the pseudoscalar sector. This amounts to a total of 4. 

• The isospin breaking parameters A^, Ey, Ay and hy which affect the non-anomalous 
HLS Lagrangian. These represent the Direct Isospin Breaking mechanism introduced in 
this paper through the BKY mechanism. 

• Some parameters floTl allowing the n° — r] — 7/ mixing. The r\ — v( mixing angle 9 P and 
the parameters named above e and e', which may account for, respectively, the 7r° — 77 and 
7T° — 7/ mixings. The last couple of parameters is not important for g — 2 estimates but 
may provide interesting physics information. One may anticipate on fit results by saying 
that the condition 9 = is well accepted by the data as in previous analyses [|38l ; as 
a matter of consequence Op can be (and will be) chosen as entirely fixed by the nonet 
symmetry breaking parameter A (see Eqs. (1261)). One will also see that the pair e and e' 
can be safely replaced by a single free parameter [61 j. Therefore, the number of really 
free parameters accounting for the n° — 77 — rf will be reduced to one. 

• Some subtraction parameters (8) involved in the mixing functions of vector mesons, in 
the p meson self-energy and in the 7 — V transition amplitudes. 

• Some more parameters (4) describing the mass and width of the narrow cu and <p mesons. 
As a detailed description of the loop corrections to their inverse propagators is of little 
importance for the present purpose , there is no need to go beyond. 

Stated otherwise, only the parameters A^, Ey, Ay and h v are new and all others have been 
already dealt with in previous releases of the present model Il39ll46ll24l . 

One may be surprised to face a so large number (~ 25) of parameters to be fitted from data. 
This only reflects that the number of physics pieces of information and of processes to account 
for is also exceptionnally large : more than 900 data points, six annihilation channels 
tt°7, 777, K + K~,K°K°, 7T + 7r~7r°), 10 radiative decay modes (VP'-f, P77 including now the 
7T° — > 77 partial width), the <\> — > ix + n~ decay and finally the dipion decay mode of the r 
lepton. All these pieces of information should get simultaneously a satisfactory description. 
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It should be stressed that the parameter space is sharply constrained, as will be confirmed 
and illustrated by the present study. One should also stress that the ix + n~ , 7r°7 and 777 cross 
sections, together with the decay modes referred to above, allow already a good determination 
of all fit parameters except for two : c\ — c 2 and A a- The former is derived from fitting the 3 
pion cross section, the second from fitting both KK annihilation channels. Actually, in order to 
accurately determine Ey, the dipion spectrum in the decay of the r lepton also plays a crucial 
role. 

This peculiarity leads us to a motivated critical analysis of the available 7r + 7r~7r°, KK and 
r data sets. As one plans to motivate a value for the hadronic contribution to g — 2, our dealing 
with the corresponding data should also be motivated. 

As far as cross sections are concerned, it is already known from our previous studies that 
the 7r + 7r~, 7r°7 and 777 annihilation cross sections are very well described within a simultaneous 
fit including also the decay data already listed. This can be seen in 113911461 : indeed Figure 2 in 
[1391 and Figures 1 and 2 in f|46l are indistinguishable from what is derived in the present study. 

10 Reanalysis of the 7t + 7t~7t° Annihilation Channel 

Taking into account the error described in Footnote |25] the analysis of the model description 
of the 7r + 7r~7T° data is worth being redone. We take profit of this case in order to exemplify 
how the dealing with data sets is done. 

The available 3-pion data sets can be gathered into 3 different groups : 

il The former data set collected by the Neutral Detector (ND) at Novosibirsk and published 
in ll9~3l : we include in this group the few data points from [1941 . These mostly cover the energy 
region in between the cu and </> peaks. 

ii/ A CMD-2 data set covering the cu region [75] together with a corresponding SND data 
sample [|9T| which actually extends up to 980 MeV. 

iii/ Several CMD-2 data sets covering the <\> region and extracted from [|8~8l |89] |90], accom- 
panied by a data set from SND [|92l starting at 970 MeV. 

The small data sample from DM1 [|95l is used for illustrative purposes and is not included 
in the fit procedures. It would not influence the fit results. 

In fit procedures, it is very hard to run minuit normally because integrating the parameter 
dependent 3-pion cross section (see Eqs. (1701 and (1711 ) renders prohibitive the execution time. 
Therefore, we still use here the iterative method described and motivated in Section 10.3 of 
061. 

The choice of the 3-pion data sets considered in the global fit was performed in [[461 re- 
lying on the data sets listed in i. Indeed, the n + n~ data used in the global fit serve to fix all 
parameters, except for the cu and <\> mass and width parameters which are derived from having 
included the 7r°7 and 777 cross sections; therefore, the ND data having a large lever arm (see 
downmost Figure|5]), they are alone able to determine accurately the value for c\ — c 2 (see third 
line in Figure |3]). 

Here one proceeds otherwise in order to learn more as each of the just above mentioned 
data set carries intrinsically a value for c\ — c 2 . Nevertheless, the group of data sets needed in 
order to fix all parameters except for c\ — c 2 has been enlarged : Beside the 7r°7 and 777 

cross sections, we have included the r decay information from ALEPH, Belle and CLEO. This 
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will be justified later on. On the other hand, one assumes C3 = C4 which is justified in Section 
US 

Fits are performed by including either the CMD-2 data sets or SND data sets, each in 
isolation. On the other hand, separate (and independent) fits are performed in either of the to 
and regions. Therefore, in these fits, the to region fits are not influenced by the <p region 
information and conversely. Moreover, CMD-2 and SND data are not influencing each other. 
The data sets associated with the so-called to and </> regions is not ours; it has been performed 
by the experimental groups who published the corresponding data sets separately. 

It should be stressed, especially in the present case, that the notion of data set covers, as 
importantly, the data points, the full error co variance matrices (i.e. including the correlations 
reflected by the non-diagonal entries), and all the additional pieces of information provided 
by the experimental groups. Among this last kind of information, the global scale uncertainty 
included in the systematics should be suitably accounted for. As far as scan data are con- 
cerned, the statistical methods we use are the standard (text-book) methods briefly reminded 
in Subsection 19.71 

The results of these fits are summarized in Figure [Hand are commented on now. As a word 
of caution, it should be noted that the experimental errors shown in these plots are the quadratic 
sum of the reported statistical and systematic errors, neglecting all correlations. As the error 
bars do not (and cannot) take into account the correlations, they should only be considered as 
a visual indication of what is going on. The real distance of data points to its best fit curve is 
instead accurately reflected by the x 2 values which, indeed, take appropriately into account all 
the reported pieces of information about the error covariance matrix. 

Top left Figured] shows the fit of only the CMD-2 data in the <p region; this provides a good 
fi@ (x 2 /npoints = 110/80) returning c\ — c 2 = 1.21 ± 0.10. Top right Figure \T\ shows the 
case for the SND data in the region in isolation; the fit is much better (% 2 /npoints = 26/33) 
but returns c\ — c 2 = 2.18 ± 0.13. These two fit values for c\ — c 2 differ by ~ lOer, clearly 
tagging an inconsistency between the CMD-2 and SND data sets in the region. 

On the other hand, one has performed likewise for the to region in isolation. One then gets 
for CMD-2 data a large x 2 1 'npoints = 26/13 with c x - c 2 = 1.29 ± 0.04 (bottom left Figure 
[T]). A closer examination of these data shows that an important part of this relatively large x 2 
is due to only the first data point which falls right on the vertical axis in this Figure. 

Instead, the SND to region data yield x 2 /npoints = 48/49 and c\ — c 2 = 1.12 ± 0.06 
(bottom right Figured). These two fit values for c\ — c 2 differ by ~ 3a; then, one may consider 
that the CMD-2 and SND data sets in the to region are in reasonable agreement with each other. 

One should note from fitting the SND to data set, the important effect of correlations : In 
the bottom right Figure [TJ the large distance of the (SND) data points to their fitting curve is 
compensated by the correlations in such a way that x 2 /npoints remains quite reasonable. The 
high level of compensation can be checked by computing the "diagonal" parJH of the x 2 which 

31 The numbers for \ 2 1 'npoints are the 3-pion sample contributions to the global \ 2 and its number of data 
points. One cannot provide the number of degrees of freedom as several hundreds of (other) data points are 
involved in each fit. 

32 Denoting by V the full covariance matrix constructed as explained in Subsection |9.7| the (full) \ 2 can be split 
up into its diagonal part X 2 dlag = E 4 ^( A i) 2 and its non-diagonal part X 2 non dia g = Eijtj V^ lA i A j> where 
Aj is the difference of the i th measurement and the corresponding value of the theoretical cross-section. 
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reflects the visual impression provided by the bottom right Figure [fl one gets Xdiag = 554 ! 

In addition, one has found instructive to plot the CMD-2 data together with the SND ones 
and the fit performed to the SND data in isolation. Thus, the bottom right Figured] illustrates 
that the correlations reported by SND allow a reasonable reconstruction of the cross section 
valid for both the SND and CMD-2 data sets. 

For information, the fit performed using only the ND dataEU yields y 2 / npoints = 25/37 
and c\ — c 2 = 1.30 ± 0.06, in good accord with the previous fit result c x — c 2 = 1.17 ± 0.07, 
derived under comparable conditions (see second data column in Table 3 of [1461 ): the difference 
between these two estimates for c\ — c 2 can be attributed to the influence of the r data samples. 

The various estimates for c\ — c 2 derived from our fits are gathered in Figure fusing obvious 
notations. Using the fit values for c\ — c 2 , as tag of consistency, this plot clearly shows that the 
region SND data set behaves differently from the other three-pion data sets. 

From these considerations, one can conclude that : 

• In the to region, there is a good agreement between CMD-2 and SND data from within 
the filter of our model. 

• In the </> region, at minimum x 2 , one can get a reasonable description of both CMD-2 
and SND data, but with much different values for the fit parameters as reflected by their 
c\ — c 2 values. 

Therefore, one observes a qualitative difference between all CMD-2 data and the SND data in 
the u) region, on the one hand, and the SND data in the <p region, on the other hand. 
One has pushed a little further the analysis by two more series of fits : 

• One has simultaneously submitted to fit the CMD-2 and SND data but only in the <p re- 
gion. One gets the result shown in Figure |2] The fit might look reasonable {x 2 /npoints = 
176/113) and returns c\ — c 2 = 1.94 ± 0.07, close to the SND value, as can be seen from 
Figure |3] 

• One has submitted separately to fit the CMD-2 data and the SND ones but simultaneously 
in the cu and <\> regions. The CMD-2 data return x 2 /npoints = 136/93 with c\ — c 2 = 
1.31±0.04, while the SND data return x 2 /npoints = 102/82 with c x -c 2 = 1.23±0.06. 
Figure @]displays the corresponding best fit curves with data superimposed^. Even if the 
X 2 /npoints and the fit value for ci — c 2 are reasonable, top right Figure @] leads us to 



avoid using the SND <\> region dat 

From this series of fit, one can conclude that it is possible to fit simultaneously the CMD-2 
and SND data in the <\> region and get a reasonable solution. However, mixing the to and <p 
regions returns, in the case of SND, an unacceptable solution, even if the x 2 /npoints may look 
reasonable. 

33 As reminded above, this data set covers the region in between the u> and cf) peaks. 

34 Here also, one may wonder that the top right Figure corresponds to a quite reasonable fit quality. We thus 
remind once more that, for all figures shown, the effects of correlated uncertainties is not - cannot be - shown. 
In the case of SND, this is larger than 5%. Along the same lines, one should mention that the errors plotted are 
always the quadratic sum of statistical and uncorrelated systematic uncertainties. 

35 In this case, the so-called "diagonal" part of the \ 2 at minimum is larger than 1 100. 
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Therefore, one is led to include in the set of data samples submitted to the global fit all 
3-pion data referred to above, except for the SND region data set. The corresponding fit has 
been performed and is shown in Figure |5]with c\ — c 2 = 1.18 ±0.03; the 3-pion data contribute 
to the global fit with x 2 / npoints = 220/179. The result shown at the last line in Figure |3] 
shows that the global fit performs, as expected, a good (fitted) average of c\ — c 2 . This also 
indicates that the data sets considered are statistically consistent with each other. 



11 Analysis of the KK Annihilation Data 

As reminded in subsection 19.51 above, several data samples are available collected by the 
CMD-2 and SND Collaborations on VEPP-2M at Novosibirsk. The CMD-2 data are extracted 
from ll88l |9~71 l98l and the corresponding SND data from [|68~1 . The quoted systematics are 
treated as a scale uncertainty and dealt with as explained in Subsection |9.71 

The published data being cross sections, the fitting function is : 



< s ) = -^T^kWk^)? with qK = -y/ s -4m? K , (K = K±,K°/K°) (78) 

for each of the 2-kaon annihilation channels; the kaon form factors F%(s) are given by Eqs. 
(|49l . Both cross sections are corrected for the intermediate photon dressing. Moreover, for the 
charged kaon channel, the additional Coulomb factor ||67l |471 Z(s), reminded in Eq. (|53l , is 
understood and is not "renormalized" as in [|68~ll69~1 with Z(mV). 



11.1 Fitting the KK Data 

In order to perform this analysis, we have done a first series of fits using separately the 
CMD-2 neutral and charged KK channels and the corresponding data from SND. In order to 
avoid peak information not following from the KK data, we have decided to remove the data 
from the 7r°7 and 777 annihilation channels from the fit procedure. However, anticipating on 
our final results, we have included the three data sets from ALEPH, Belle and CLEO referred 
to above. 

Therefore, the additional data sample is composed of all e + e~ — > n + n + data, all — > 
tx^tx^v data and 18 partial width decays (all VPj and P77 modes and also the three leptonic 
decays (p/cu/(p) — > e + e~ modes). None among these pieces of information has any direct 
influence on the description of e + e~ — > KK, even through the <\> mass and width parameters 
which are, thus, solely determined by the KK data. 

The results are shown in Figure |6] left side for the K°K° data and right side for K + K~ . 
One observes a good description of the K°K° data for each of the samples provided by the 
SND or CMD-2 Collaborations. The picture is quite different for the K + K~ data; the CMD-2 
data sample is well fitted, while the SND sample is poorly fitted. Additional information for 
these peculiar fits is displayed in the first two lines of Tabled One can see there, that the value 
for x 2 /N associated with the K°K° data are the same for both data samples, while they differ 
significantly for the corresponding K + K~ data samples. 

Fitting simultaneously both CMD-2 and SND K Q ~K° data samples only, returns the same 
X 2 information, illustrating that the corresponding data samples are perfectly consistent with 
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each other. Simultaneous fits of all KK data confirm this property (see third line in Table |2). 
Interestingly, the % 2 's at best fit in the third and fourth lines practically coincide with the sum 
of the corresponding information in the first two lines of the same Table. This illustrates that 
the so-called additional data set sharply constrain the KK cross sections. Moreover, in view of 
the fit results for CMD-2 data, one can consider that the constraints are well fulfilled by data, 
giving a strong support to our modelling. 

The ratio of cross sections a(e + e~ — > K°7( )/o-(e + e~ — > K + K~) is observed to provide 
a valuable piece of information, as it allows to magnify the effects mentioned just above. This 
is shown in Figure |7J where the data for this ratio are plotted normalized to the ratio of cross 
sections as coming out from our fits. The data ratio plotted in the top Figure [TJis derived from 
the information given in [[69] and one can estimate its uncertainty to ~ 2.3 4- 2.4 %. 

The CMD-2 data points normalized to the fit expectations bin per bin is perfectly consistent 
with 1 over the whole s region covered by the <p resonance. The dotted lines in top Figure [7J 
represent the experimental scale uncertainty and do not take into account the uncertainties on 
the fitting functions. This also illustrates that our modified Breit-Wigner lineshape is very well 
accepted by the data. 

In contrast, the SND data exhibit a behavior reasonably well averaged by the fit function 
ratio; however, it does not look consistent with flatness - at least as well as for CMD-2 data. 

It follows from these considerations that the largest self-consistent data set for the KK 
channel is made by merging all CMD-2 data and the K !^ data provided by SND (see last 
line in Table [2]). 

As a matter of information, beside getting an appropriate description of both e + e~ — > 
K°K^ and e + e~ — > K + K~ cross sections, it is worth remarking that the radiative partial 
widths included in the fitted data set are also well accounted for. For instance, including also 
the e + e~ — > (n° /rj)^ cross sections in the fitted data set, the remaining set of 10 radiative 
decays yields a quite remarkable x 2 /n = 6.5/10, with estimated r(7r° — > 77), T(r] — > 77) and 
r(V — >■ 77) at respectively 0.27a, 1.77<r and 0.23a from their accepted values [10]. As the 
corresponding couplings are strongly affected - especially g n o 71 - by (see Eqs. (1581). we 
consider that physics validates our model. 



11.2 The HLS Solution of ^ KK Puzzle 

The partial width decays — > KK are defined by : 



Sir, 



T{<p^KK) = ^-\g 4>K - K Y 
Therefore, one has : 



(q. 



K 



Am 2 K ) 



(79) 



TU ->• K + K~ 



T(d> -> K°K 



Br(0 -> K + K~ 
Br(0 -> K°K C 



R 



9(j>K+K- 



Z(ml)~RZ(ml)(l-2A A ) (80) 



where R = 1.528 originates from the ratio of momenta and the Coulomb factor computed at 
the 4> peak is Z(m^) = 1.049. The ratio of couplings has been given in Eq. (l54l . Therefore, 
using A a from the last line in Table |2] one gets : 
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Br(0 -> K+iT-) 

— ^- zzn^ = 1.41 ±0.03 (81) 

Br(0 -> K°k) 

The same ratio can be computed from information given by CMD-2 in a recent paper f[69l 
and amounts^ to 1.47 ± 0.04. The difference between the CMD2 estimate and ours amounts 
to about 2a. Our final result, obtained by using the largest possible ensemble of data sets, 
provides A^ = (6.34 ± 0.70) 10 -2 and then the ratio of branching ratios becomes 1.40 ± 0.02. 

Therefore, the HLS model, equipped with the (BKY) direct isospin symmetry breaking 
mechanism, provides a solution to the long-standing puzzle concerning the <p ~* KK decays 
as thoroughly analyzed in ||47l and more recently discussed in [|68l . In our approach, the 
mechanism responsible for this is, in fine, the kaon field renormalization which should be 
performed within the HLS model once isospin symmetry breaking is performed a la BKY- 
BOC. Indeed, as the neutral and charged kaon field renormalization factors play in opposite 
directions (see Eq. (fl~8l ). they pile up in the ratio. 

The relatively large value found for A^ indicates that several sources contributes to the 
BKY breaking of isospin symmetry. The contribution to A^ due to the light quark mass mass 
difference ll58l (~ 1%) is certainly not the single source and others - like electromagnetic 
corrections - are certainly absorbed within the numerical value for A^. Moreover, it is also 
likely that different corrections at the VK + K~ and VK !^ vertices may influence the fit value 
for A a- Being global, the BKY breaking mechanism cannot allow to disentangle the various 
contributions to A^ which share a common order of magnitude (each at the percent level). The 
situation is quite different from the breaking of SU(3) symmetry which is widely dominant 
numerically and can motivatedly be compared to ChPT expectations 11381 . 



12 Analysis of the r Decay Data 

Using F£(s), the pion form factor in the decay of the r lepton (see Subsection 16.21) . the 
partial width of the two-pion decay is given by Eq. (|73T) . On the other hand, the quantity which 
encompasses the full experimental information in this field is Eq.(l72l) : 

_ 1 dT n7T {s) 1 dN(s) 

H(s) - -—j—S EW G EM (s) - B^--^- 

as, indeed, the lineshape and the absolute magnitude of each experimental spectrum are merged 
together. The full width r r is taken from the RPP ifTOll . The last two factors in the mid- 
dle expression above account for isospin symmetry breaking effects specific of the r de- 
cay : Sew — 1-0235 for short range corrections 111031 . G E m(s) for long range corrections 

In former studies, it was shown that the lineshape alone was perfectly consistent with an- 
nihilation data ll39l |24| . However, if one also takes into account the absolute magnitude - 
represented by the branching ratio B W7r in the formula reminded just above - the agreement 
is poor. In order to reach a satisfactory description of the data, Ref. Il24l introduced a mass 

36 The uncertainty might be somewhat overestimated, as one has assumed independent the errors for Br(</> — > 
K+K~) and Br(0 K°K°). 
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difference Sm 2 and a coupling difference Sg between the neutral and charged p mesons, which 
underlays all reported stand-alone fits to r spectra [13] . However, additional scale factors 
were needed and their fitted values were found consistent with the reported scale uncertainties 

However, the present study, as reflected by Table [TJ above, has clearly demonstrated that 
isospin breaking of the HLS model does not necessarily result in non-vanishing Sm 2 and Sg at 
leading ordeQ As emphasized above, the BKY-BOC breaking scheme instead leads to a dif- 
ference between the universal vector coupling (g) as it comes in the anomalous sector and in the 
non-anomalous sector of the HLS Lagrangian (g(l + £y)). We prove, here, that this provides 
a much better account of all data than only assuming some mass and width differences supple- 
mented with some residual rescaling. Stated otherwise, it is because Direct Isospin Breaking 
acts differently in the anomalous and non-anomalous sectors that the model yields an almost 
perfect description for all data, without any need for some additional rescaling. In this mecha- 
nism, the single sensible difference between the pion form factor in e + e~ annihilations and in 
r decays resides in the difference between the transition amplitudes 7 — V and W — V. 

Figure [8] shows the global fit result for the function H(s) defined just above together with 
the data points from ALEPH 00), Belle ED and CLEO 02) Collaborations B The inset 
magnifies the p peak region. One can clearly conclude to a nice agreement between model and 
data, all along the fitted region - from threshold to 1 GeV. The corresponding pion form factor 
in e + e~ annihilations coming out of the global fit is represented in Figure |9j These two Figures 
illustrate that the simultaneous description of e + e~ and r data allowed by the model is, indeed, 
as successfull in both sectors. 

Figure [10] shows in two different manners the r residual behavior. Top Figure [TOl displays 
the usual residuals for the function H(s), while downmost Figure ITOl represents (H mode i(s) — 
Hdata(s)) I H mode i (s). These can be compared with respectively Figure 3 and Figure 4 from Il24l 
where the (Sm 2 , Sg) parametrization of isospin breaking was used. The comparison clearly 
indicates that the present model better performs for all r data sets and, especially, for the 
ALEPH flU data. 

In order to allow for a deeper comparison with the previous release ll24l of the present 
model, we reproduce in Table [3] (first data column) the fit results reported in [|24l together with 
our new fit results under various conditions. 

The second data column in Table [3] is derived excluding the KK data sets in order to be 
as close as possible to E4l . One observes, for almost all data sets, better fit results than in the 
former release of our model [24 j . There is no effect in introducing the 3-pion data set from 
SND I19T1 (covering the to region) as the xln/dof = 1.11 is unchanged. It is also worth noting 
that the partial width for 77 — > 77 is found at 0.43a from its accepted value [fTOl : the distance is 
O.llcr for rf 77 and 0.47a only for the newly introduced n° — > 77 decay mode. 

One may conclude therefrom that the HLS model, equipped with the mixing schemes pro- 
vided by loops and by the direct isospin breaking procedure, provides a fully satisfactory solu- 
tion to the e + e~ — r puzzle, both in magnitude and in shape. The relatively poorer fit quality 
for the BELLE data might be related with the absolute scale issue revealed by the stand-alone 

37 0ur present results as well as formers [24] tend to indicate that an electromagnetic correction to the p mass 
does not give a significant effect (see Footnote|6]l. 

38 When dealing with r plots, the error bars represent the diagonal errors, i.e. no account of bin-to-bin correla- 
tions is attempted. 
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fiH] provided by BELLE ||4~T1 . Therefore, one can confirm that : 

• The main drawback of the breaking model in [|24l was a too tight correlation between the 
universal coupling in anomalous and in non-anomalous processes. This has been cured 
by defining the Direct Isospin Breaking mechanism substantiated by a highly significant 
value for £ v = (3.74 ± 0.42)%. 

• The breaking model in lfi~3l may account insufficiently for the difference between the 
p° — 7 and — W ± transition amplitudes. 

Therefore, the reported discrepancies between the pion form factor in e + e _ annihilations 
and in r decays can always be attributed to an incomplete treatment of isospin symmetry break- 
ing. For information, Figure [TT] displays the ratio of the transition amplitudes / P7 and f p w as 
coming from the global fit and already given in Table [1] : 

/; h v A v a(s) V2/3(s) 

7? + ^ _ + ~ + ^ _ ^ 

We have found appropriate to provide in the third data column of Table[3]the results of the fit 
obtained keeping the KK data sets, while excluding all the data sets. The fourth data 

column reports on the fit quality reached using the full data set we considered safe. This means 
all data sets discussed above, except for two SND data sets : The e + e _ — > 3n data set collected 
above 970 MeV ||92l and the e + e~ — > K + K~ data set. These have been shown to provide either 
an unacceptable behavior for the fit solution [92], or a poor x 2 Il68l . In this configuration, one 
fits 906 data points (including the 10 individual decay modes) corresponding to 881 degrees 
of freedom. The global fit probability is highly favorable (71%). This configuration will be 
referred to in the following as "Solution A" or "Configuration A". 

In this Solution A, one observes some tension between the KK and data groups. 

Indeed, comparing its content with the second data column, one observes that the data 
group yields a x 2 increased by 30 units. Instead, comparing Solution A with the third data 
column in Tabled one does not observe any significant degradation the fit quality of the KK 
data group : The x 2 f° r the K°K° data group is improved by 2 units, while the x 2 f° r K + K~ 
data group is worsened by 6 units. 

As this 30 unit increase of the (7r + 7r~7r°) x 2 ma y l°°k abnormal, we have tried tracking its 
origin. This issue is clearly related with having introduced the KK data which influence the 
model description of the <fi region. Therefore, we have redone fits excluding all the 7r + 7r~7r data 
data sets covering the 4> region. One obviously remarks a significant effect; this configuration 
will be named hereafter "Solution B" or "Configuration B". 

In the following, any differential effect between what has been named Solutions A and B is 
examined carefully and commented. 

39 The fit published by BELLE reveals a very significant improvement if the absolute normalization of their 
spectrum is left free; instead of returning an absolute scale of 1, the best fit exhibits a significant ~ 2% shift. 
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13 Structure of the u — > tttt Coupling 



As noted in Subsection I6.3L the coupling u — > nn in the upgraded broken HLS model is 
given by : 

<W = y [(1 - h v )A v - a(s)} . (82) 

This expression exhibits two contributions of different origin. The first part is a constant term 
generated by the Direct Isospin Breaking procedure defined at the beginning of this paper, 
the second is generated by the kaon loop mixing procedure already defined in [|39l l46l and 
reminded above. This structure resembles that given in [|481|49|| . It is interesting to examine the 
behavior of the ratio : 

^ = [(1 - h v ) A v - a(s)\ 

as a function of y/s. It is given in Figure[[2l where the vertical line figures the to mass location. 
Of course, the effective part of this function is determined by the tu Breit-Wigner distribution 
and is concentrated within a few tens of MeV's apart from the to peak position. 

From the best fit discussed in the above Section (see the second data column in Table [3]), 
one gets the central values for the fit parameters and their error covariance matrix. These have 
been used to generate g W7r7r by Monte Carlo methods. Computed with using the RPP ifTOl mass 
for the to meson, this gives3 : 

' g uw7T = (-0.071 ± 0.003) + z'(0.150 ± 0.002) 

< A y = (-5.22 ±0.75) 10" 2 , h v = 1.690 ± 0.107 (83) 

a = 2.288 ±0.006 , g = 5.556 ± 0.014 , £ y = (3.74 ± 0.50) 10~ 2 

The observed useful correlations are < ST,ySAv >= —0.056, < ST,y5hv >= 0.028 and 
< 5h v 5A v >= 0.232. 

In order to stay consistent with [|4~8l l49l definitions, one can consider that g 1 = ag(l + 
£y)/2 and g^ nn = ag(l — hv)Av/2 are the couplings of the ideal fields, defined as such before 
applying the loop mixing. Therefore, the quantity G : 

G = 4^ = (1 - h v )A v (l - E v ) (84) 

9 pTTTT 

should be close to the parameter carrying the same name in [|49Tl . One finds G = (3.47 ± 
0.64) 10~ 2 to be compared with the two estimates of the same parameter given in [|49l : 
G = (7.3 ± 3.2) 10" 2 when relying on the data from [E3 and G = (4.4 ± 0.4) 10~ 2 when 
using, instead, the [176117711 data. 

Referring to Eqs. (|28l . one can conclude^ that there is much more isospin inside the 
physical p than isospin 1 inside the physical tu. In this case, one also gets for the direct term 

40 The quoted uncertainties for Ay, hy , a and XV are the improved uncertainties returned by the routine MINOS 
of the MINUIT package 11071 . 

41 The isospin component inside the physical p meson is given by ZiyAy, while the isospin 1 part inside the 
ijj is given by (1 — hy ) Ay 
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ag(l — hv)Ay/2 = —0.332 ± 0.024. Comparing this number with 7le(g U7T7r ), it is clear that 
ag(l — h v )A v /2 and TZe{a{m^)) compensate to a large extent, in such a way that g W7r7r is 



14 The 7T° — t) — rj Mixing Properties 

The mixing of pseudoscalar neutral mesons has been addressed in Section 0] especially 
in Subsections 14.31 and 14.41 The present Section is devoted to examining how the upgraded 
breaking scheme developed in this paper performs compared to the results previously derived 
in this field [38 1. In order to perform this study, we let free the pseudoscalar mixing angle 
Op, which mostly determines the relationship between the physical rj and rf fields and their 
underlying octet and singlet components rf R and rf R . The parameters e and e' which account 
for, respectively, the n° — rj and n° — r/' mixing are also let free. 

As shown in [38] and revisited in Subsection l4.3l above. the ChPT mixing angles [f50l[5TTl 0$ 
and # 8 can be expressed in terms of the nonet symmetry breaking parameter A (or, better, using 
instead v defined in Eq.([22l)). z A the SU(3) breaking parameter of the Lagrangian Ca and the 
singlet-octet mixing angle P . Therefore, they can be estimated from fitting the data already 
defined. 



The mixing angles #0 and 0$ have been recently introduced with the 2-angle description of 
the ij/rf mixing [[5011511 . The broken HLS model provides expressions for these in terms of 
the singlet-octet mixing angle P and of the breaking parameters za and A (see Eq. (|26l) and 
also [[381). 

Therefore, using the fit results (parameter central values and their error covariance matrix) 
one can reconstruct the values for and 8 . Having left free Op, one obtains the results shown 
in the first data column of Tabled Therefore, as in former studies, one observes that is small 
and its distance to zero is only 2.8<r; this should be compared with the estimate O = —4° given 
with no quoted uncertainty in iBTll . The value for 8 is numerically as expected from other 
kinds of data iBTTl . The 'tHooft parameter ll35l A is found of the order 10 %, twice smaller than 
in [|38l where an approximate treatment of nonet symmetry breaking was used. Finally, the 
singlet-octet mixing angle P is still found twice smaller than 8 , as in the former study 11381 . 

As the distance to zero of is 2.8<r, the non-identically vanishing of O is on the border 
of statistical significance. Therefore, imposing the condition O = is worth being considered; 
this turns out to algebraically relate 0p to z a and A by tan 0p = tan£> (see Eqs. ([26l) ). Per- 
forming such a fit returns the results shown in the second data column of Table [4] with a quite 
comparable probability. 

It is interesting to observe that the value for S is nearly unchanged and that the value for 
A is affected below the 10~ 4 level only. One also observes that the value for P generated by 
the appropriate Eq. (|26T ) is found in agreement with its fitted value (when this parameter is 

42 In traditional fits with the Orsay phase parametrization of the u contribution to the pion form factor, this 
property is reflected by a value for this phase close to tt/2 




14.1 The Mixing Angles Q , 9 8 And 6 P 
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left free). We conclude therefrom that assuming 9 = does not degrade the fit quality and is 
consistent with data. 

One should also note that the nonet symmetry breaking parameter A = 8.5% has a statistical 
significance of 2Aa. Performing an approximate nonet symmetry breaking f|38~1 , the value for 
A was overestimated by a factor of 2. 



14.2 The 7T° — t) and tt° — r/ Mixing Properties 

These mixing properties are reflected by the parameters named respectively e and e' as 
displayed in Eqs. (|23l) . Comparing analogous fits performed by letting free and unconstrained 
9p, e and e', we did not find sensitively different results than those obtained by imposing the 
constraint on 9 P resulting from the condition 9 = 0. Therefore, from now on, all presented fit 
results will refer to this configuration. One should note that the numerical results given in the 
above Sections have also been derived under these conditions. 

The global fit returns e = (4.89 ± 0.44) 10" 2 and e' = (1.68 ± 0.44) 10~ 2 , reflecting that 
the 7T° — r] mixing is certainly much more important than the n° — r{ mixing phenomenon. With 
the concern of reducing the number of free parameters, we have also assumed ItoTj : 

\/2 cos 9p — sin 9p 

e — e cos9p—^= 

V2 cos 9 P + sin 9 P 

(85) 



v2 cos 9 P + sin 9 P 

= — zeosinflp— p 

V2 cos9p — •m\9p 



with 9 P still determined by the constraint 9 = 0. This reduces the number of free parameters 
by one more unit. The fit returns eo = (3.16 ± 0.23) 10~ 2 with an unchanged probability; this 
corresponds to values for e and e' very close (2a each) from the corresponding fitted values, 
while the global fit probability is unchanged. The partial widths for the three decays P — > 77 
are all well accounted for : 1.64er (77), O.ller (77') and 0.06cr (tt°). Additional fit detail can be 
found in Table |3] 

The question of whether the present eo can be identified with the variable carrying the same 
name in [61 ] is uncleaS Indeed, an important part of isospin symmetry breaking effects are 
already included in the definition of the renormalized PS fields (see Eqs. (TT9b and (12TT) ) which 
undergo the rotation defined by Eqs. d23T > . Therefore, our e, e' and e carry only a part of 
the isospin breaking effects, while another part (governed by A^) has been propagated to all 
sectors of the effective Lagrangian. 



15 The Values of the FKTUY Parameters 

Our global fit modelling is in position to provide the most accurate information concerning 
the parameters c 3 , C4 and c\ — C2 defining the scales of the various FKTUY anomalous pieces 
[1281 of the HLS Lagrangian. 

43 The quantity named eo in [61] is related with R = (m s — m)/(m<j — m u ) by eo = V^3/ (4-R). For instance, 
[108| gives R = 37.2 ± 4.1, while II 1091 relying on QCD sum rules proposes R = 33 ± 6. These provide 
respectively e = (1.16 ± 0.13)% and e = (1.31 ± 0.24)%, which have little to do with our fit result. 
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In order to get the most accurate results, we have explored the parameter behavior and 
found that the least correlated combinations are c 4 + c 3 , c 4 — c 3 and c\ — c 2 . Running under the 
configuration A defined above, one gets : 

- ^ C i+ C3 - 0.962 ±0.016, c = = (-3.98ii;l) 10 -2 , c x - c 2 = 1.208iffi 



2 ' 2 

with g = 5.541 ± 0.016, while configuration B leads to : 



(86) 



c+ ee = 0.978 ± 0.020 , c_ ee = (-6.75™) 10- 2 , d - <* = 1.123j£ffi 

(87) 

with g = 5.530 ± 0.015. The correlation coefficients are similar in both cases : < [5c + ] [5(ci — 
c 2 )] >~ -0.20, < [5c + ][5c_] >~ -0.10 and < [<Jc_][«J(ci - c 2 )] >~ 0.80. Therefore, our 
global fit yields quite consistent numerical values wathever the configuration^ for the FKTUY 
parameters. 

These values can be compared with existing estimates. Using the 7r°77* form factor, [1231 
yields c + = 1.06 ± 0.13, while the partial width u — > 7r°7 provides c + = 0.99 ± 0.1— when 
using g = 5.80 ± 0.91. Our own estimates are consistent with these with, however, (MINOS) 
uncertainties five times more precise. 

A rather unprecise value for the ratio c = c_/c + has also been derived 11231 relying on the 
decay to — > ir°fi + n~, c = 0.42 ± 0.56, consistent with our results but still much less precise. 

From our results, which happen to be the most precise in this field, one may conclude that 
data only favor a partial fulfilling of the VMD assumptions f|23ll , in the sense that c 3 — c 4 = 
is in agreement with data at the 2a level, while C\ — c 2 + c 4 = 4/3 is badly violated. This can 
be rephrased as follows : the VMD assumptions [|23l are experimentally fulfilled in the triangle 
anomaly sector and strongly violated in the box anomaly sector. This confirms the previous 
parent analysis 11461 and former studies on the box anomaly in the 77/77' — > 7r + 7r~7 decays [|37ll . 

In order to go beyond, better data on the annihilation channels involving anomalous cou- 
plings ([ir°/r)}j, 7r + 7T~7r°) are needed; including new processes like the 77/77' —> 7r + 7r~7 decay 
spectra or information on the annihilation channels may also help as their dependence 

upon c 3 — C4 or c\ — c 2 is more important than in the previous channels. 

It thus follows from the present analysis that assuming c 3 = c 4 is justified. In this case, one 
obtains the following results : 

c 4 = c 3 = 0.950 ± 0.014 , ci - c 2 = 1.194 ± 0.060 , g = 5.556 ± 0.014 , (88) 
for Configuration A and : 

c 4 = c 3 = 0.951 ± 0.016 , ci -c 2 = 1.169 ±0.060 , g = 5.553 ± 0.012 (89) 
for Configuration B. 

In both cases, the correlation coefficient is < [<5c 3 ][5(cx — c 2 )] >~ —0.20. Therefore, the 
condition c 4 = c 3 drastically reduces the correlation among the surviving FKTUY parameters. 



44 Running our code excluding the KK data (see second data column in Tabled yields c + = 0.967 ± 0.021, 
c_ = (— 5.18±1;|§) 10~ 2 and a - c 2 = 1.0T4±g;g||) with g = 5.530 ± 0.015. This configuration pushes the 
significance for a non-zero c_ at the ~ 1.7a level. 
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Moreover, the fit quality is not significantly changed while assuming 04 = 03. Indeed, config- 
uration A yields x 2 /dof = 858.08/882 (71.2 % probability) instead of X*/dof = 853.98/881 
(73.7 % probability) and configuration B x 2 /dof = 728.38/882 (97.0 % probability) instead 
of x 2 /dof = 722.05/881 (97.9 % probability), where the difference mostly affects the set of 
partial widths which is always well fitted. Therefore, the improvement obtained with the up- 
graded breaking model is not due to releasing the condition c 4 = c 3 . From these considerings, 
it is justified to impose c 3 = c 4 for the rest of this study. 

16 Hadronic Contributions to g — 2 

In |[24l. one analyzed in full detail the hadronic contribution to g — 2 of most of the data 
sets used in the present study. The framework was the previous release of the present model 
studied in detail in ll46l l24l . Within this framework, only the simultaneous account of both 
annihilation channels to KK was missing. On the other hand, one might find unsatisfactory 
that some global rescaling of experimental r dipion spectra was still playing an important role, 
even if this rescaling was in accord with expectations. These two issues motivated the present 
study. 

As shown above, the upgraded model allows by itself a satisfactory account of all consid- 
ered spectra simultaneously. It is therefore worth reexamining within our upgraded framework, 
how the hadronic contribution to g — 2 is estimated and how this estimate evolves depending 
on the various kinds of data groups considered. 

16.1 The tt+tt Contribution to g - 2 : VMD Estimates 

The most important hadronic contribution to g— 2 is the n + ix~ channel. Several experiments 
[175117611781 and some analyses |[T3ll24l[TT0ll give the n + ix~ contribution to g — 2 provided by the 
energy region [0.630 — 0.958] GeV. Therefore, it is worth considering the information provided 
by this reference region; this allows to substantiate the improvement which can be expected 
from VMD-like models. Indeed, several kinds of information are worth considering : 

• While unifying the description of e + e~ annihilation and r decays, one expects an in- 
creased precision on the anomalous magnetic moment of the muon a^(nix). 

• While having a framework which encompasses most of the physics up to the <p region, 
the stability and the robustness of the a^nn) estimates can be examined. The relative 
statistical consistency of the various data groups is also an issue which can be addressed, 
relying on their behavior under global fits. 

Table |5] displays our estimate for the nir contribution to = (g — 2)/2 provided by the 
reference energy range under various fit configurations. In each case, the fitted (central) pa- 
rameter values and their error covariance matrix are used in order to sample several thousand 
parameter vectors, assuming a n-dimensional gaussian error distribution. Each vector of sam- 
pled parameter values is, then, used to compute a^(nn). The corresponding distribution of the 
<V(7T7r)'s is then fitted to a Gaussian function. The results displayed in Table [5] are the central 
values and the standard deviations of this distribution which intrinsically takes into account 
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the correlations among the fitted parameters. Unless otherwise stated, the FSR correction is 
included in all reported contributions of the n + ix~ channel to a M . 

Beside the experimental spectra, there is always a set of partial width decays submitted to 
fit. These have been defined in Subsection 19.61 In the results reported below, one should keep 
in mind that the accepted values [10] for the (p/w/0) — > (vr°/?7)7 and — > e + e~ partial 
widths are included in the set of partial widths submitted to the fit as long as the experimental 
spectra for the e + e~ — > (n° /rj)~j annihilation channels are not used. As emphasized in 11461 . 
this hides some model dependence which might be somewhat conflicting with our own model. 
This explains why one should prefer any configuration where the e + e~ — > (w° ' /rf)^ data are 
submitted to the global fit. 

Also, when the data for the two annihilation channels e + e~ — > KK are not considered in 
the fit, one chooses to fix eo = = 0, as we have no real sensitivity to them. Likewise, 
ci — c 2 is absent from fits as long as the e + e~ — > data are not considered. Finally, 

the parameters fixing the mass and width of the <p meson are left free only when the fitted data 
allow to constrain them. 

In the first line of Table [5J one finds the value for a^ix-ix) derived by submitting to fit the 
scanned data for the annihilation process e + e~ — > n + n~ - together with the full set of partial 
width decays. This result compares well with the value derived using the previous release of 
our broken HLS mode£f[ as can be seen by comparing with the relevant piece of information 
reported in Table 4 of [1241 . 

As there is no longer any mismatch between e + e~ and r data, both in magnitude and in 
lineshape (see Section [12]), it is legitimate to merge them. This merging provides the new and 
important result given in the second line of Table [5] One clearly observes that the merged 
nn data give a result perfectly consistent with the e + e~ — > n + ix~ data alone with a quite 
nice probability. The central value for a^ixix) is nearly unchanged and the uncertainty slighty 
improved. 

This is, of course, the main effect of having upgraded our symmetry breaking procedure of 
the HLS Lagrangian. In this new framework, there is no need for an auxiliary rescaling [1241 of 
the r spectra and the net result is a perfect consistency of the e + e~ — > ix + n~ data with/without 
the r data considered as constraints. This statement can be substantiated by comparing this 
result with those reported in the entry "NSK+ A B C" of EU (a^irir) = (364.48 ±1.34) 1CT 10 ) 
which exhibited a shift of about 5 1CT 10 produced by the three r data sets, a ~ 3.6a effect. 

The third line in Tabled displays the effect of replacing the (p/co/(p) — > (n° /rj)~j and the 
(u)/(f)) e + e~ partial widths by the cross sections for e + e~ — > (ir°/r))j. The central value 
for a^-Kix) is practically unchanged, while its standard deviation is increased by ~ 9%. The 
following line in Table [5] displays the effect of including the full e + e~ — > data group 

already defined. As in 11241 . one observes a perfect consistency of the results for a^nn). In 
total, the standard deviation is slightly reduced (a(a^(inr)) ~ 1.5 1CT 10 ). At this point one 
may conclude that the central value is marginally modified by fully including the (p/w/0) 
(n°/r])^ and e + e~ — > data groups within the fit procedure. The variations of the 

uncertainty returned by the fits might rather reveal statistical fluctuations. 

45 Even if expected, this proves that the effects produced by having introduced Sy do not modify the fit descrip- 
tion of the e+e~ — > tt + tt~ data. 
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The last line in Table [5] displays the effect of including the two e + e~ — > KK cross sec- 
tions into the fitted data set. One observes some effect, as a^(nn) undergoes a 1.9 1CT 10 shift 
upwards while the fit probability remains quite good. This fit configuration - referred to as 
Solution/Configuration A - encompasses the largest set of data samples considered safe. This 
turns out to consider that the 30 unit increase of the x 2 associated with the data group, 

even if large, is not abnormal (see the fourth data column in Table [3]). 

The result shown in the last line of Table [5] may reveal some tension among the data set 
groups. In order to explore this issue, one has redone fits excluding the n + n~n data group, and 
examined the effects of using the selected e + e~ — > K°K° and e + e~ — > K + K~ data, either 
separately or together. The corresponding results are displayed in Table [6J Comparing the 
statistical information here with those in the last line in Table \5\ renders somewhat suspicious 
the quoted 30 unit increase of xl+n-jr - 

A final piece of information is provided by performing the fit using the data group 

data amputated from the data points collected in the region above 1 GeV (therefore, excluding 
the 4> region). This fit configuration has already been referred to as Solution/Configuration B. 
The reason which motivates this removal is that the 7r + 7r~7r data before introducing the KK 
data is only constrained in the <fi region by the relatively unprecise data on the 7r°7 and r^y 
channels. One then obtains : 

f Solution B : a^vrvr) = (362.44 ± 1.49) 1(T 10 x 2 /dof = 722.05/801 Prob. = 97.9% 

[ Solution A : a^(vrvr) = (362.19 ± 1.44) 10~ 10 x 2 /dof = 854.00/881 Prob. = 73.7% 

(90) 

where the result for Solution A is reminded. 

These differences indicate that all physics channels covering the region are worth to be 
reconsidered, as already argued from discussing the fit results in TableO Indeed, the difference 
in fit quality between Configurations A and B reveals some tension between the KK data and 
the 7r + 7r~7T° data collected in the <\> region. Fortunately, the physics in the <\> region is still 
accessible at VEPP-2M. It seems also in the realm of the KLOE detector, as this turns out to 
run DA$NE within a ± 20 MeV interval apart from the 4> mass peak value. 

16.2 The 7r + 7r Contribution to g — 2 : Comparison with Data 

An interesting piece of information comes from comparing our (VMD) estimates derived 
from global fitting with the corresponding estimates provided by the various experimental 
groups. 

Table [7J displays the published experimental results concerning the contribution of the 
0.630 — 0.958 GeV/c region to a^nn). We first list the three important results from CMD-2 
and SND; as we also use the data sets from OLYA and CMD [80], we also give at the line 
flagged by "OLD" our average using these data sets together with those from NA7 HI 111 . TOF 
Eia, M2N [El, DMl EQ, all collected before those from Il75ll76ll77ll78l. 

The third data column provides, first, our average derived using the data sets from [17311761 
[TTJITHl and, next, also those including the older data sets referred to just above. Our results are 
directly comparable with these as we do not yet use ISR data. 
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Both solution A and solution B results favorably compare with the scan (mi) data averaging 
as the uncertainty is reduced by a factor close to 2. 

The following lines of Table [7] display, for information, the experimental results derived 
from the data sets collected using the ISR method and the global average of the ISR and scan 
data. 

One should stress that our results for a^irn), derived excluding the ISR data, provide 
information already comparable in precision to those obtained using them. This motivates to 
examine the ISR data in view of including them into the fit procedure. 

One may also compare our estimates with the weighted average of the r data HOl |4~T1 |42| 
which gives 10 10 a /J (7T7r) = 365.21 ±2. 67 exp in the reference region, including FSR corrections; 
applying the p — 7 corrections proposed in |[T6l . this becomes 10 10 a M (7T7r) = 361.66 ± 2.67 exp 
and provides 10 10 a /i (7T7r) = 361.15=1=1.76^ when averaged with the e + e~ data. This indicates 
that examining the idea proposed in lfT6l in a wider context is an interesting issue. Indeed, this 
could lead to another successful VMD-like model and, therefore, may contribute to a motivated 
evaluation of the model dependence of a M estimates. 

As a summary, one may conclude that our global model provides a good determination of 
the contribution to a^(irn) from the invariant mass region 0.630—0.958 GeV/c. The accuracy of 
our VMD estimates is found much improved compared to direct averaging of the experimental 
data and their central values are found consistent within uncertainties. By including ISR data 
at a later stage, the precision of the result might be further increased. 

16.3 Hadronic Contribution to g — 2 

In Table [8l one displays the contribution of each of the examined channels to a M from their 
respective thresholds up to 1.05 GeV/c, i.e. slightly above the </> peak. 

The first two data columns show the results corresponding to the so-called configura- 
tions/solutions A and B. These have been derived by fitting the data sets referred to in the 
preceding Sections and the motivation to consider both solutions valid can be emphasized from 
Tabled 

The last two data columns exhibit the averages of experimental data for each of the mea- 
sured channels submitted to the global fit. These differ by excluding (third data column) or 
including (fourth data column) in the averaging the ISR data sets collected by KLOE fTTl [T9ll 
and BaBar [|T8l for the n + ix~ final state. As we have excluded for now the ISR data from our 
analysis, the gain due to the global fit can be directly inferred by comparing with the third data 
column; nevertheless, it is interesting to compare the accuracy of solutions A and B to the 
averages derived using the high statistics ISR data. 

As expected, the improvement generated by the global fit affects all the channels considered 
and is always a factor of 2 or more (see the ix + n~n Q channel) better than the average of the same 
data. The first line even shows that our accuracy is comparable - actually slightly better - than 
the average derived using the ISR data. 

It is interesting to note that the sum of all contributions for solution B is in accordance 
with the result expected from the standard sum as reported in the third (or fourth) data column. 
Solution A, instead, gives a smaller sum than the experimental average of the same data; the 
distance is 2.97 10~ 10 , i.e. ~ 1.6cr^ eor . or ~ 0.7a exp . 
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It is interesting to examine the individual channel contributions. Those from the 7r°7 and 777 
channels, as calculated from data, rely on pretty poor statistics and generally cover restricted 
energy ranges ll82l l83l l84l l87l l86l l85l (see Subsection 19.21 ); instead, our model results are 
estimated (significantly) larger and cover precisely the full energy range from thresholds to 
1.05 GeV. This especially concerns the region in between the u and <p peaks. 

Our model estimates for the 7r + 7r~7r° and K + K~ channels are found smaller than the ex- 
perimental averages at the 1 or 2 a exp levels, while the K°K contribution corresponds to the 
experimental expectation. This confirms the need for a better experimental knowledge of all 
annihilation channels in the <\> region. 

The first data line in Table |9]reports the results derived from fits with our global model. The 
second line ("missing channels") provides the experimental averaged contribution to from 
the channels unaccounted for within our model (the Air, 5n, 6ir, r]nn and tun final states). This 
has been computed using the trapezoidal integration rule. As the corresponding data are sparse 
below 1 .05 GeV, this estimate might have to be improved. 

The line "Total Model" provides the estimate of the full hadronic vacuum polarization 
(HVP), merging our model results with the additional listed contributions. 

The corresponding experimental average taking into account all available ISR data sets 
EH QH ED has been estimated 0U to a M (e+e~) = (690.75 ± A.72 tot ) 1(T 10 , including the 
contributions above 5.2 GeV calculated using perturbative QCD. For comparison, the corre- 
sponding total average provided by HI 101 is a^(e + e~) = (695.5±4.0 ea ;p±0.7Q C <£)) 1CT 10 (not 
accounting for the recent KLOE data set lfl~9l ); accounting for all the available ISR data sets, 
lfl4l yields as experimental average a fJ ,(e + e~) = (692.3 ± 4.2 toi ) 1CT 10 . 

In order to illustrate the impact of r data, we present separately the fit results derived when 
including or when excluding the r data sets from the fitted data sets, keeping for the rest the 
configurations leading to solutions A and B as previously defined. 

Including r data sets results in an increased value of the hadronic VP by ~ 3 1CT 10 . This 
will be commented on below. One also remarks that our uncertainties are comparable to the 
experimental one, even if our estimates are penalized by having - provisionally - discarded the 
ISR data. Our estimates also compare favorably with the revised estimate excluding all ISR 
data given by El : a^e+e") = (690.9 ± 5.2 exp+rad ± 0.7 QCD ) W~ 10 . 

16.4 The Anomalous Magnetic Moment of the Muon a M 

Table [lOl displays our final results concerning a M . We still report on the results derived in 
the fit configurations A and B, using or not the r data in the fit procedure. The leading-order 
(LO) hadronic VP discussed in the previous Subsection is reminded in the first line. In order 
to yield our estimate of a M under the various quoted configurations, one should add the effect 
of higher-order hadronic loops taken from I1T61 , the light-by-light contribution [5|; we took 
the latest estimate of the pure QED contribution^ and the electroweak (EW) contribution 
is taken from [4]. Summing up all these, one obtains the values given as "Total Theor." which 

46 The recent [114] value a^QED] = 11658471.8096(0.0044) displayed in Table [TO] should be updated to 
a^lQED] = 11658471.8960 (in units of 10 -10 ). In order to compare with already published results we prefer 
keeping the former value for our estimates of the HVP and of g — 2. 
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should be compared with the average [1 J of the different measurements for a M , recently updated 

m. 

The difference between our theoretical estimates and the experimental average [2] is finally 
given together with their respective statistical significance. The significance of this difference 
varies between 4.07cx (solution B including r's) to 4.65cr (solution A excluding r's). The 
difference between including r's and excluding them is a ~ OAa effect. [D] provides an 
estimate excluding the KLOE data lfTTll - and the more recent ISR data sets not available at 
that time - reaching a difference with the BNL average [2] of (30.1 ± 8.6) 10~ 10 , a 3.5a 
significance. Our least significant estimate (solution B including r's) is, instead, 4.07cr. 

Figure [T3] displays our results together with the most recently published estimates. On top 
of the Figure, one finds the estimates using or not the r data provided in [fl~4|. The following 
entry is the estimate given in [fi~6l which combines e + e~ and r data (after correcting for the 
p° — 7 mixing). The last entry H115I is derived including the ISR data (HLMNT1 1); this is the 
latest result using the final KLOE El and BaBar HH data. 

We have also displayed the latest result [13] derived excluding ISR data which directly 
compares to ours. This indicates that the improvement provided by the global fit method cor- 
responds to increase the discrepancy of the BNL measurement [2] with the Standard model 
prediction by ~ 0.6 -r- 0.8a. Therefore, the discrepancy starts reaching an interesting signifi- 
cance. 

16.5 Influence of Data Set Choices on the Estimate for 

In order to derive our estimates for a M , we have defined a paradigm, unusual in this field. 
Indeed, one usually performs the average using all data sets contributing to a given final state 
in isolation; the prescription used is the S-factor technics of the Particle Data Group. However, 
this supposes the simultaneous handling of statistical and systematic uncertainties. The most 
common way of performing this handling is to use as weights the quadratic sum of statististical 
and systematic uncertainties |]9). 

In our approach, especially in this paper, the underlying paradigm is different and can be 
formulated in the following way : 

• All different channels are correlated by their underlying common physics and an Ef- 
fective Lagrangian approach is presently the best tool to deal with the non-perturbative 
QCD regime. 

• All data sets, covering or not the same physics channel are considered by taking into 
account the peculiarities of their uncertainties as reported by the experimental groups. 
There is, in principle, no real difficulty in order to deal with statistical uncertainties. It 
is commonly assumed that uncorrected systematics and statistical uncertainties could 
be added in quadrature and we followed this rule. Other systematics involving bin-to- 
bin or experiment — to-experiment correlations should be treated as such; the method is 
standard q| and has been sketched in Subsection 19.71 

47 In the scan experiments we deal with in the present paper, all reported correlated systematics can be consid- 
ered as global scale uncertainties for which the standard method applies. For ISR experiments lfT7l[T9l [T8l. the 
situation is different as several independent sources of systematics are defined which, additionally, vary all along 
the spectra. The standard method can be extended to this case [46]; however, it should better be reformulated in 
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• The Lagrangian model should allow for a good description of a large number of data sets 
in as many different physics channels as possible. The goodness of the global fit should 
be accompanied by a good description of each group of data sets - ideally each data set. 
As tag for this property, we choosed the x 2 / n points value for each data set group; this tag 
should not too much exceed 1 . Referring to our case, the 7r + 7r~, 7r°7, 777 physics channel 
data and the reported partial width decays already represent an acceptably good start- 
ing point, allowing a critical examination of the data associated with further additional 
channels. 

• Including a new data set, or a new group of data sets, should not result in a significant 
degradation of the already accounted for data sets. This should be observed at the global 
level and at the local levels (i.e. for each group). Following from the analyses in Sections 
[TOl and [TH peculiarities of their fit behavior led us to discard from our global fit the 
K + K~ data set and one of the 7r + 7r~7r° data sets provided by SND. This turns out to 
require that the (large) set of data samples considered be statistically self-consistent : 
Only 2 data sets out of 45 did not pass this consistency criterium. 

At this point, given the (broken) Lagrangian one uses, the selection criteria are only the 
global fit quality and the "local" (data set specific) fit properties reflected by the various x 2 /n points 
values, discarding any possible consequence for the value for a M . With Solutions A and B, one 
has also avoided any kind of data set reweighting by discarding the two data sets exhibiting 
some faulty behavior compared to the rest. 

Nevertheless, it is a simple exercise to switch on the two discarded SND data sets within our 
fitting code. For information, this leads to Aa M = (a At ) eip — (a^) th = (34.00 ± 8.21) 10~ 10 , 
a 4.14a effect. However, this is associated with an exceptionally poor global fit probability 
(1.75%) and to xl+ n -^/n points = 331/212 = 1.56 and xl+K- /joints = 93/62 = 1.50. 
Interestingly, and somewhat unexpectedly, the x 2 / n points for the other data sets are practically 
unchanged compared to Table |3] except for the decay data set account which is sharply de- 
graded : Xdecays/ n points = 20.5/10 ~ 2. This may reflect that our broken HLS model is so 
sharply constrained that poor data sets are mostly reflected by poor global fit probabilities. 

A tag value of x 2 l n po%nts = 1.3, as yielded for the chosen final state data, is on the 

border of what could look reasonable to us (see third data column in Table [3]). Nevertheless, 
compared with x 2 / n points = 11 (see second data column in Table [3]), it looks acceptable; 
however, this corresponds to an increase by 30 units of the absolute magnitude of ^+ , 
when introducing the selected kaon data. One may, indeed, consider that this indicates some 
tension within the 4> region data calling for a closer experimental examination which can be 
performed at the existing facilities covering the <p region. 

Awaiting for better data in the region, we have been left with two challenging solutions : 
Solution A which uses all the data sets we have considered as secure, and solution B obtained 
by removing all 7r + 7r~7r° data sets above the KK threshold. 

a way which avoids introducing as many scale factors to be fitted as sources of different systematics. Indeed, this 
may produce fit instabilities and, on the other hand, one has to deal with correlations between physics parameters 
and these scale factors which may be uneasy to handle. 
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16.6 The Differential Effect of the Various r Data Samples 



In view of the discussions above, we have chosen to display all our final results for a^, in 
the fit configurations corresponding to solutions A and B. On the other hand, as can be read off 
Table |3l at the fit properties level, one can consider that the so-called e + e~ — r puzzle is over. 

However, one still observes a (2-j-2.5) 10~ 10 increase of the returned values for a p produced 
by the r data. As stated already above, the r data are essential in order to return a reasonably 
precise value for our fit parametei@ Sy. Therefore, the shift attributable to the r data can be 
considered as a normal consequence when fitting a model with a more constraining set of data 
samples. 

Nevertheless, Table |3] indicates that the x 2 / n points are sensitively different for ALEPH (~ 
0.43), CLEO (~ 1.26) and BELLE0(~ 1.77). This difference of fit quality leads us to examine 
the effects of removing the CLEO data sample and/or the BELLE data sample for our fitted data 
set. 

When keeping only the ALEPH data sample, we get Aa M = 38.47 ± 8.22 (a 4.68cr signifi- 
cance) and Aa M = 36.81 ± 8.90 (a 4.13cr significance) for respectively solutions A and B. As 
can be seen from Table [101 these strikingly resemble the corresponding values for Aa^ derived 
when keeping only e + e~ data in our fit procedure (i.e. excluding all r data). In these peculiar 
configurations, the ALEPH data fit quality which was already very good (% 2 /n points ~ 16/37), 
becomes impressively better (x 2 /n points ~ 4/37). 

Going a step further, we have examined the effect of considering only ALEPH and CLEO 
data. In this case, our fit returns Aa M = 36.02 ± 8.22 (4.38<r significance) and Aa^ = 34.74 ± 
8.26 (4.21cr significance) for respectively solutions A and B. One can check with Table [101 that 
these values become closer to their partners when fitting excluding r samples. 

Therefore, using only the r data samples from ALEPH HOI and/or CLEO 11421 returns 
values for Aa^ consistent well within errors with those derived using only e + e~ data. The 
slightly different behavior of BELLE data may be related with the normalization issue sketched 
in footnote |49] 

16.7 On the Significance of the HLS Value for Aa M 

In view of the considerations developed in the two preceding Subsections, one can certainly 
consider that the most conservative estimates for Aa M are those derived while including r 
data as they are reported by ALEPH, BELLE and CLEO. This corresponds to the information 
provided in the first two data columns of Table [TO] 

This means that the disagreement between the BNL measurement [J2] and the Standard 
model prediction for Aa M lays in beween 4.07 and 4.33 a. Moreover, from our analysis of the 
differential effects of the various available r data samples, one may consider these bounds as 
conservative and that the significances in the right part of Table [101 cannot be discarded. 

48 The numerical accuracy of the scan e + e~ data alone does not permit a precise determination of Sy which 
is returned by MINUIT with large errors. 

49 Leaving free the absolute normalization of their dipion spectrum improves the stand-alone fit of the BELLE 
Collaboration BP from 80/52 to 65/51. This corresponds to a best normalization of 1.02 ± 0.01. Such a re- 
normalization of their absolute scale has some influence on the value for a M . One should remind that we do not 
have any longer fitted rescaling factors in our fitting functions. 
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In view of this, in the perspective of taking into account relatively poor data setgroup, one 
has rerun our code in order to get the solution when weighting the contributions ofj : 

• all 7T + 7r~7r° data in our global sample by 179/232.41, 

• the BELLE data sample by 19/32.31, 

• the CLEO data sample by 29/36.48, 

in the global x 2 while leaving the other weights (all equal 1) unchanged. This turns out to 
rescale globally the uncertainties associated with the corresponding data sets by the inverse of 
these weights, assuming that their relatively poor quality is only due to an overall underestimate 
of the uncertainties by a factor of respectively 1.14 (tt+tt^tt ), 1.30 (BELLE) and 1.12 (CLEO). 
This may look as a way to infer some sort of S-factors inside the global fit procedure. 

This reweighting procedure^] provides as total hadronic VP contribution to a M (686.32 ± 
4.60) 1(T 10 and Aa M = (34.93 ± 8.23) 1(T 10 , a 4.25a significance. 

Going a step further, another check may look appropriate. As the contributions of the 
7r + 7T~7r°, BELLE and CLEO data to the total x 2 have been weighted in order to reduce their in- 
fluence, one can do alike with those groups of data which exhibit too favorable individual x 2 's. 
Still referring to fitting with configuration A, this turns out to weight the "Old Timelike" data 
by 82/56.61, the tt°7 data group by 86/68.37, the 777 data group by 182/123.3L the ALEPH 
data by 37/15.92 while keeping unit weights for the "New Timelike" and both KK data groups. 
This leads to an hadronic VP of (685.00 ±4.58) 10- 10 andto Aa^ = (36.25 ±8.21) 10" 10 cor- 
responding to a a 4.41a discrepancy. This is almost identical to the value found with Solution 
B, excluding r's, as can be seen from Table [TOl 

Therefore, these exercises enforce our conclusion that the most conservative value for Aa^ 
exhibits a discrepancy of 4.07a and values as large as ~ (4.30 + 4.50)cr are not unlikely. 



17 Conclusion and Perspectives 

Several aspects should be emphasized. They can be grouped into two items : Low energy 
hadronic physics description and g — 2 related topics. 

Concerning the first item, the present study indicates that the HLS model suitably broken 
is able to encompass most low energy physics in an energy range extending up to the_d> me- 
son mass. More precisely, among the non-baryonic possible final states, one coverso most 
channels with multiplicity n < 4. 

More precisely, equipped with the so-called upgraded direct symmetry breaking - in the 
u, d and s sectors - and including the mixing of neutral vector mesons produced at one-loop, 

50 The weights used in this Subsection refer to partial x 2 's obtained by fitting under Configuration A with 
assuming C3 =04; it is the reason why they slightly differ from the corresponding numbers given in Table [3] 

51 We have also made a fit leaving free scale factors affecting the covariance matrices of the 3-pion data as a 
whole, of the BELLE and CLEO data. The hadronic VP we get is (686.73 ± 4.49) 1(T 10 , quite similar to this 
value. 

52 Among these, only the process e + e~ — > rjirir has not been examined; however, the good description of the 
77/77' — > 7T7T7 decays reported in [46 1 indicates that it could be successfully considered. On the other hand, the 
e + e~ — > ujir annihilation is too much influenced by high mass vector resonances [29 , 30] to be accounted for by 
the standard HLS model. 
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the HLS model accounts quite satisfactorily for all the examined physics pieces of information. 
This covers the 6 annihilation channels having significant cross sections up to the meson 
mass and a few more spectra like the dipion spectrum in the r decay and, also, an additional 
list of partial width decays. Previous studies [|37l [461 have also shown that the dipion spectra 
in the 77/7/ — > mvy decays fall inside the scope of the HLS model. 

It is an attractive feature of this framework to exhibit a parent character between the long 
reported issues represented by the e + e~ — r and the </> — > KK puzzles : Indeed, it is the 
same breaking mechanism implemented in the Ca and in the Cy pieces of the HLS Lagrangian 
which provides a solution to both. It permits - together with the s-dependent vector meson 
mixing - to finalize the consistency of the e + e~ and r physics and to reproduce the branching 
fraction ratio <p K + K~ /(f) — > K°K° . This is materialized by a satisfactory simultaneous fit 
of both e + e~ — > KK cross sections and of the pion form factor in both e + e~ annihilation and 
r decay. 

The upgraded model thus provides a tool allowing a simultaneous treatment of a large 
number of experimental spectra. It also permits a critical analysis of the fit behavior of any 
data set in consistency with the others. Then, one is in position to discard motivatedly some 
data samples which do not behave satisfactorily within a global fit procedure and could then 
put some shadow on derived numerical results. We have shown that such data samples are 
only few : 2 out of the 45 considered spectra. It should be stressed that discarded data sets are 
always identified because of their full redundancy with some other data sets, which are found 
to behave normally within the global model; stated otherwise, this removal is not expected to 
produce a bias and, a contrario, any effect resulting of keeping them is suspicious. 

The model provides a tool which has the virtue of exhibiting the physics relationship be- 
tween the various physics channels. Within the global fit procedure involving the data on each 
channel, the model parameters yield a better accuracy which propagates to all the reconstructed 
pieces of information, especially the photon hadronic vacuum polarization and, thus, improves 
significantly g — 2 estimates. 

Indeed, we have shown that the various components of the HVP yield central values in 
accordance with expectations and an uncertainty improved by a factor of 2 quite uniformly 
within the fit range. This has been shown for the 7r + 7r~, 7r°7, 777, 7r + 7r~7r°, K + K~ and K°K~° 
channel contributions up to 1.05 GeV. Up to this energy, these channels represent altogether 
more than 80% of the hadronic VP and one of the two dominant sources of uncertaint\o 




In order to figure out the gain in terms of statistics, one can make the following state- 
ment : considering globally the existing data sets is equivalent to having x4 more statistics 
simultaneously in each of the considered channels without any increase of the systematics. 
Therefore, considering additionally the high statistics ISR data leaves some room for improved 
estimates of the HVP, provided the dealing with systematics can be reasonably well performed. 
One should nevertheless stress that the global method we advocate, used with only the standard 
scan data samples provides already as good results as all scan and ISR data using the standard 
numerical integration of the experimental cross sections. 

One may also try to figure out the improvement expected from including the high statistic 
ISR data samples [fT71[P9l[T8l within the fit procedure. Being optimistic, one may think that the 

53 The other dominant error comes from the hadronic VP between 1 .05 and 2 GeV. 
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uncertainty on the HVP contribution up to 1.05 GeV could be divided by 2, from ~ 2 x 10 
(see Table [8]) to ~ 1 x 1CT 10 . Let us also assume that the ISR data samples will not rise 
unsolvable bias problems. Taking into account the rest of the HVP, which carry an uncertainty 
of ~ 4 x 10~ 10 (see Table |9]), the uncertainty on the full HVP would decrease from ~ 4.60 x 
10~ 10 (see Table |9]) to ~ 4.25 x 10~ 10 . Using the information collected in Table [KB the total 
uncertainty on would decrease from ~ 5.30 x 10~ 10 to ~ 5.00 x 10~ 10 and the uncertainty 
on Aa M would decrease from ~ 8.20 x 10 -10 to ~ 8.00 x 10~ 10 . This may look a marginal 
improvement; the reason for this is the large value for the systematics generated by hadronic 
HVP in the region 1.05 -j- 3.10 GeV (see Table©, which thus becomes a prominent issue for 
future significant improvement^ 

However, this is not the end of the story. In the course of the paper, and this is well expressed 
by Tables l9l and [TOl we saw that below 1.05 GeV systematics may produce significant shifts of 
the central values for the HVP and thus for a M . This was observed, for instance, in the A and 
B configurations, where the shift for the HVP - and for a M - amounts to ~ 2.00 x 10~ 10 (see 
also Subsection I16.6I ). Because of this, there is still valuable experimental work to do also in 
the sub-GeV domain to decrease and/or better understand systematic errors. More precisely, 
a better experimental knowledge of all channels in the mass region - 0.95 4- 1.05 GeV - 
may result in improving quite significantly our estimate on g — 2 and in resolving some of the 
ambiguities discussed in the main text. As stated above, the information in this mass region 
has an important influence down to the threshold regions. This is certainly within the scope of 
existing machines and detectors^!. 

What are the prospects for the future? 

A new muon g—2 experiment at Fermilab is expected to come into operation in 5 years from 
now. The accuracy is expected to improve to 0.14 ppm from its current 0.54 ppm. This also 
requires a factor 4 improvement of the hadronic vacuum polarization. As demonstrated by our 
analysis, it is possible to improve the low energy part up to and including the by a systematic 
application of effective field theory methods in form of a resonance Lagrangian approach. 
However, as mentioned above, the main effort will be required in the range above the (f> up 
to about 3 GeV. In this range, major progress is expected from CMD3 and SND at VEPP 2000 
at Novosibirsk, from BESIII at Beijing, as well as from exploiting additional yet unanalyzed 
ISR data from BaBar and Belle. Within the 5 years available until a new experimental result for 
a M will be realized, lattice QCD is expected to be able to produce results which are competitive 
with standard evaluations based on data. This also would provide important cross checks for 
the present results and, more generally, for the effective Lagrangian approach. 

For now, one can conclude that the paradigm represented by a global model which encom- 
passes the largest possible set of data indeed results in a highly significant improvement of the 
photon HVP uncertainty and of the uncertainty on g — 2. As the global model allows to detect 
problematic data sets susceptible of generating biases, it must be accompanied by the most 

54 Actually, even if the uncertainty on the HVP contribution coming from the energy region up to 1.05 GeV 
vanishes, this would not entail a significant improvement of the global uncertainty for a M ! Stated otherwise, 
reducing the HVP error in the region from threshold to 1.05 GeV from ~ 4 x 10~ 10 to ~ 2 x 10~ 10 has much 
more dramatic effects than reducing it from ~ 2 x 10~ 10 to ~ 1 x 10~ 10 . This is a pure algebraic effect following 
from having to perform quadratic sums for final uncertainties. 

55 One may remark that scan data for the e + e~ — > tt + it~ cross section in the region are still not available. 
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accurate possible treatment of the reported experimental systematics. 

Taking into account the ambiguities generated by a limited number of data sets, the most 
conservative estimate for the hadronic vacuum polarization leads to a significance for a non- 
zero Aa M of 4.1a. Solving these ambiguities discussed in the main text may result in a signifi- 
cant increase of this conservative bound. 
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Appendices 



A The Full HLS Non-Anomalous Lagrangian before Loop 
Mixing 

The non-anomalous Lagrangian of the Hidden Local Symmetry Model can be written : 

C-HLS = (£ A + C V ) = C-VMD + C T (91) 

in order to split it up into convenient pieces. Removing the pseudoscalar field kinetic energy 
term, which is canonical, one has : 
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(92) 

in terms of the first step renormalized vector fieldso. The pseudoscalar fields shown here are 
renormalized (it is the origin of the za and terms). Of course, we have only kept the lowest 
order symmetry breaking contributions. 

Some parameters have been introduced in Eq. (|92l for convenience; these are (m 2 = 

ag 2 f 2 n ) : 



2 2 2 H i ^ 1 2 2 

m p o = m u = m [1 + ^vj , m <i) = m zy 



fp-y = agfn 



1 + Sy + h V - 



«5/ 2 



l + £y + 3(1 -/iy)Ay] , 



(93) 



On the other hand, using : 



m 2 ± = m 2 [1 + Ey] , f pW = agfl [1 + £ 



(94) 



56 In order to avoid heavy notations, the subscript iZj, which actually affects each of the vector fields in Eq. 
has been removed. 
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one has at lowest order in the breaking parameters : 
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(95) 

where one has limited oneself to write down only the terms relevant for our purpose. The 
(classical) photon and W mass terms [|23ll34l are not considered and have been given only for 
completeness. However, it is worth remarking that the photon mass term does not prevent the 
photon pole to reside at s = as required [52], at leading order. 

Our breaking scheme generates new couplings for the charged p mesons : 
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(96) 

because of the field redefinition given by Eqs. (fT9l , (|22l) and (1231) . Therefore, the broken 
HLS model predicts decay modes r — > ■n(r}/rf)v of small intensity absent from the original 
Lagrangian. 



B Elements of the SM 2 Matrix 

The perturbation SM 2 to the full mass matrix M 2 is defined in Eq. (1321 . Keeping only the 
leading terms in isospin breaking parameters, its matrix elements are : 
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(97) 



The functions ei(s) and £2(5) and the constant g pK K have been already defined in the main text 
by Eqs. (1331) . We have also defined : 



IpKK - 4 



(98) 



C Lagrangian Pieces with Renormalized Vector Fields 

Coupling to a pion pair comes from the two Lagrangian pieces^ : 
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(99) 
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which exhibit the couplings to a pion pair depending on mixing angles. 

Similarly, the Lagrangian pieces relevant for couplings to K + K~ are given by : 
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and by : 
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for K°K couplings. Setting b = a(z v — l)/6 and p = z v V^, the s-dependent loop transition 
functions IL/ 7 are : 
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where 



es(s) = e 2 (s) + ei(s) and e D (s) = e 2 (s) - e^s) 



(103) 



57 Throughout this Section, one takes profit of introducing irrelevant second-order terms in breaking parameters 
in order to write down expressions in the most concise way 
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The expressions in Eqs. (11021) are very close to their partner in [|39l or [|4~6l . as only first-order 
perturbation terms are meaningful. 



D The Anomalous Lagrangian Pieces 

The full Anomalous Lagrangian can be written : 

^-anomalous = £>VVP + C-AVP + £>AAP + C-VPPP + C-APPP 



(104) 



where A denotes the electromagnetic field. It incorporates the Wess-Zumino-Witten terms and 
the FKTUY Lagrangian f28ll . The Lagrangian pieces occuring in Eq. (11041) are0 [23] : 
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(105) 
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o , -[1 - -( Cl - c 2 + c 4 )]e^^Tr[g^P9 Q P^P] 



where the q are parameters not fixed by the model. iV c is the number of colors fixed to 3. The 
V and P field matrices are the bare ones. 



E The Vr x P^ Coupling Constants 

In order to express the V^Pj couplings, it is appropriate to define the angle 5p = Op — 6$ 

(tan e = 1/V2) : 

sin#p = —j=(cos 5p + v^sin^p) 

f (106) 

cosOp = —j={y2 cos 5p — sin^p) 
v3 

and some parameter expressions which reflect the various ways, nonet symmetry breaking in 
the PS sector occurs : 
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For clarity, the new constant parameters are denoted exactly as they are defined in 
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where v is the nonet symmetry breaking parameter defined in Eq. (|22l i. Finally, we also have 
defined G = — eg(c 3 + c 4 )/(87r 2 / 7r ). The p Q Rl P^i coupling constants are : 
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In the cup 1 P'j sector, one has : 
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and, finally, the (pp^P'j sector provides much simpler expressions : 
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Finally, the sector is described by : 
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where z T is another breaking parameter [|39l l46l not discussed here. 
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F The VPPP Coupling Constants 

The VPPP coupling constants in the P n ~ix + (P = ir°, 77, r/') have been defined for the 
Ri renormalized fields in Eq. ||62T ). With an obvious naming, they are obtained by multiplying 
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by D = —3g(ci — c 2 — c 3 ) / (An 2 fl), which depends on the FKTUY parameters c x — c 2 and c 3 
not constrained by the model. Only the leading correction terms have been retained. 



G The VPj Couplings for Renormalized Vector Fields 

Let us define the quantities : 



egc 3 , 
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for each Vr 1 = pp 4 , w^, $p x and P = tt , ??> The fi , y J?i P o7 can be found in AppendixE 

in Eqs. (11081) . (11091) . (11101) . The functions Hy° occuring in Eq. (|66T > provide the couplings of 
the physical vector fields to a photon and a neutral meson. They are given by : 



(114) 



These definitions help in writing the cross sections in a way similar to those in 11461 . When 
expanded, the Hy° functions may contain contributions of order greater than 1 in some of the 
breaking parameters. These higher-order contributions are irrelevant and can be dropped out. 
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H The Functions NAs) in e + e — > n°n + n Annihilations 



The amplitude for the transition 7* — > is much simply expressed in terms of the 

following complex functions : 
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2e 1 (cos6 P - V2sm6 P ) + 2e 2 (\/2 cos 6 P + sin6» P ) - A A \/3 1 



6^/3 



D p o(s) 



^ + Ms) - (1 - hy)A V ] %M + 7( S )^T 

Do{s) D^s) 



D u {s) 
1 



1 1 

+ 7= — ; r + 



D p0 {s+-) D p +(s -) D p -(s 0+ ) 



N 3 (s) = [a(a+_) - (1 - V)Av] 



1 



iV 4 (s) 
N 6 (s) 



ei(cos^p — \/2sin0p) + e2(v / 2cos#p + sin#p) 
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^ - (1 - h v )A v 



1 



-D p o(s+_ 



7(«+-) -7(g) 
2L» p o( S+ _) 

2ei(cos 6>p - y/2 sin 6>p) + 2e 2 (V2cos6>p + sin6>p) - A A \/3 1 



6\/3 



D p o(s+_) 



(115) 

s + _, s - an d s 0+ are the invariant mass squared of the corresponding pion pairs from the final 
state, s is the off-shell photon invariant mass squared. All other parameters and functions have 
been defined in the body of the text. 

The connection with the Kuraev-Siligadze (x, y) parametrization [|73l is defined by(m = 

m w o, m-x = m n ±) : 

s + _ = s(2x + 2y — 1) + ml 



s+o 



s{l-2y)+ml 



(116) 



s_o = s(l — 2x) + ml 

The integration limits can be found in [73]; they are also reminded in 11461 . The Kuraev- 
Siligadze kernel is : 

2 

I (H7) 



G(x, y) = A(x 2 -)(y 2 - — ) - \ 1 - 2x - 2y + 2xy + 

s s V 



2m? — m, 



7T U 
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Xl/N 


Xl/N 


A A (%) 


Fit Prob (%) 


K°K° + K+K- (SND stand-alone) 


60.10/60 


56.54/26 


8.54 ± 1.93 


33.7 


K°K° + K+K- (CMD-2 stand-alone) 


59.30/59 


29.00/36 


5.98 ±0.86 


85.8 


K°K° (SND & CMD-2) 


115.68/119 




5.51 ±3.21 


81.8 


K°K° + K+K- (SND & CMD-2) 


119.83/119 


88.09/62 


6.29 ±0.80 


40.4 


K°K° (SND & CMD-2) 
+ K+K- (CMD-2) 


118.54/119 


29.27/36 


6.09 ±0.79 


80.8 



Table 2: Fit quality of the K+K~ and if if data. Beside the additional data sample (see 
text), each line in the first column tells which KK data samples have been included in the fit 
procedure, xl is the \ 2 value for K°K° data, xl is the corresponding information for K+K~ 
data. The iV's are the respective numbers of data points. The last data column provides the 
global fit probability for each case. 
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X 2 /N 


(5m 2 , Sg, c 3 = c 4 ) 




Statistical Information 






m 


excl. KK 


eXCl. 7T + 7T 7T° 


A 


B 


Decays 


16.20/9 


5.53/10 


6.13/10 


11.36/10 


5.94/10 


New Timelike n + ix~ 


126.47/127 


119.73/127 


130.33/127 


127.50/127 


129.65/127 


Old Timelike n + ix~ 


60.45/82 


51.64/82 


56.36/82 


56.09/82 


56.60/82 


7T°7 


66.07/86 


66.84/86 


61.19/86 


67.21/86 


66.93/86 


777 


135.78/182 


128.89/182 


122.64/182 


122.62/182 


121.37/182 




139.44/126 


200.92/179 




230.98/179 


105.91/99 


K+K- 


— 


— 


29.93/36 


35.16/36 


29.85/36 


K°K° 






120.07/119 


117.94/119 


119.99/119 


ALEPH 


36.51/(37+1) 


21.25/37 


15.92/37 


16.80/37 


16.16/37 


Belle 


28.29/(19+1) 


27.02/19 


34.19/19 


32.22/19 


33.62/19 


CLEO 


39.46/(29+1) 


35.12/29 


35.86/29 


36.09/29 


36.03/29 


X 2 /dof 


648.68/680 


656.93/726 


612.63/703 


853.98/881 


722.05/801 


Global Fit Probability 


80.1% 


96.8% 


99.4% 


73.7% 


97.9% 



Table 3: Comparison of the fit qualities between the fit results of the model as it was in (|24| 
(second data column) and as it is now (third data column). KK data were not submitted to fit 
in [|24l . The '+1' added to the number of data points for r data stands for the experimentally 
given r.m.s. affecting the (fitted) global scale. The 3-pion data set information is displayed 
boldface in order to show the difference in the fit data set : In the second data column, the 
3-pion data set from SND [|9T| has been (newly) introduced and in the last data column only 
the 3-pion data sets collected below the region are considered. 
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General Fit 


Constrained Fit 




-1.11° ±0.39° 







-23.88° ± 0.34° 


-23.82° ± 0.34° 


9p 


-12.66° ±0.35° 


-12.91° ±0.18° 


A 


(8.52 ±3.55) 1(T 2 


(8.52 ±3.55) 10~ 2 



Table 4: Some parameter values derived when leaving free 6p and A (first data column) or 
when relating them by imposing 9 = to the fit (second data column). 



Data Set 


Fit Solution 


Statistical Information 






X 2 /dof 


Probability 


e + e~ — t> 7r + 7r~ 


360.00 ± 1.64 


177.38/208 


93.3% 


+ [r] data (ABC) 


359.8 ±1.47 


262.94/293 


89.6% 


++ (e + e _ -> /rf\~f) 


360.09 ± 1.60 


436.94/549 


99.9% 


++ (e+e~ ->■ 7r + 7T _ 7r°) 


360.91 ± 1.45 


661.22/727 


96.1% 


++ (e+e" ->■ M) 


362.79 ± 1.43 


858.08/882 


71.2% 



Table 5: The contribution to 10 10 a M (7r7r) from the invariant mass region 0.630 — 0.958 GeV/c. 
The first line provides the fit results using all the e + e~ — > it + it~ annihilation data set group. 
The next line uses the previous data group and the three r spectra. By "++" at any given 
line, we always mean all data sets belonging to the groups referred to in the preceding lines, 
plus the data set group indicated at this line. FSR corrections are taken into account. An 
appropriate set of radiative decays is always understood. The last line refer to what has been 
named Solution/Configuration A. 
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Pit Solution 


y 2 /dof 

A / U<J1 


Probability 


OTllv P~^~ P — ^ fC 
will y c c r J \. .i v 


360 79 + 1 4Q 


474 69/585 




only e+e- ->■ 7T 7Z° 


362.83 ± 1.47 


580.78/668 


99.34% 


both e+e" -> 


362.81 ± 1.47 


613.29/704 


99.40% 



Table 6: The contribution to 10 10 a M (7r7r) from the invariant mass region 0.630 — 0.958 GeV/c 
using KK data sets under various conditions. All ix + ix~ix° data have been excluded from fit. 
FSR corrections have been performed. 



Data Set 


Experimental Result 


Average 


Fit solution 


CMD-2 (1995)11751 


362.1 ± (2.4) sta * ± (2.2) syst 






CMD-2 (1998)117511771 


361.5 ±(1.7) stat ±(2.9) S2/st 


SND (1998) 11781 


361.0 ±{1.2) stat ±{4.7) syst 


Average 




361.26 ±(2.66) tot 


A : 362.79 ± 1.43 tot 


OLD 


354.1 ±(3.3) stat ±(8.1) syst 




Average (excl. ISR) 




360.65 ± (2.55) iot 


Fit Solution A 
Fit Solution B 






B : 363.16 ± 1.47 tot 


KLOE-2008 Ca 


356.7 ±(0.4) stat ±(3.1) S2/st 






KLOE-2010 030 


353.3 ± (0.6) sta4 ± (3.2) syst 


BaBaR fISllllOl 


365.2 ±(1.9) stat ±(1.9) syst 


Total Average 




360.53 ± (1.44)** 



Table 7: The various published estimates of the contribution to 10 w a fl (7iir) from the invariant 
mass region 0.630 — 0.958 GeV/c. The quoted averages always refer to all experimental 
results displayed in the preceding lines. The line "OLD" information refers to our average 
performed using the data sets collected before those of CMD-2 and SND (see text). Our fit 
solutions A and B are derived using the r spectra from f|40ll42ll4TI . KLOE-2010 estimate for 
a^nir) is ours, as the experimental spectrum stops slightly below y/s = 0.958 GeV [fT9l . 
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Process 


Solution B 


Solution A 


Data (excl. ISR) 


Data (incl. ISR) 


7T + 7T~ 


498.54 ± 1.97 


497.98 ± 1.76 


498.53 ± 3.73 


497.72 ±2.12 


7T°7 


4.64 ±0.04 


4.28 ± 0.04 


3.35±0.11 tot 


VI 


0.65 ±0.01 


0.67 ±0.01 


o.48 ± om tot 


rfl 


0.01 ±0.00 


0.01 ±0.00 




7T + 7r~7r 


42.03 ±0.60 


40.88 ±0.52 


43.24 ±1.47^ 


K+K- 


16.87 ±0.20 


16.93 ±0.18 


17.88 ±0.54 tot 


K°K° 


12.02 ±0.09 


12.07 ±0.08 


12.31 ±0.33 tot 


Total Up to 1.05 GeV 


574.76 ±2.10 


572.82 ± 1.90 


575.79 ±4.06 tot 


574.98 ± 2.66 tot 



Table 8: Contributions to 10 from thresholds up to 1.05 GeV/c The experimental er- 
rors merge the reported statistical and systematic uncertainties in quadrature. FSR effects 
(3.43 10 -10 ) have been included into the 7r + 7r~ contribution. The first two data columns dis- 
play our fit results and the last two data columns report the direct numerical integration of the 
relevant data. 
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Final State 


Range (GeV) 


Contribution (inch r) 


Contribution (excl. r) 






Solution A 


Solution B 


Solution A 


Solution B 


e + e~ — > hadrons 


threshold ->■ 1.05 


572.82[1.90] 


574.76[2.10] 


569.86[2.15] 


571.40[2.27] 


missing channels 


threshold ->■ 1.05 


1.55(0.40)(0.40)[0.57] 


J/rj) 




8.51(0.40)(0.38)[0.55] 


T 




0.10(0.00)(0.10)[0.10] 


hadronic 


(1.05,2.00) 


60.76(0.22)(3.93)[3.94] 


hadronic 


(2.00,3.10) 


21.63(0.12)(0.92)[0.93] 


hadronic 


(3.10,3.60) 


3.77(0.03)(0.10)[0.10] 


hadronic 


(3.60, 5.20) 


7.64(0.04)(0.05)[0.06] 


pQCD 


(5.20, 9.46) 


6.19(0.00)(0.00)[0.00] 


hadronic 


(9.46, 13.00) 


1.28(0.01)(0.07)[0.07] 


pQCD 


(13.00,oo) 


1.53(0.00)(0.00)[0.00] 


Total 


1.05 ->■ oo 
+ missing channels 


112.96 ±4.13** 


Total Model 


threshold — > oo 


685.78 ±4.55 


687.72 ±4.63 


682.82 ±4.66 


684.36 ±4.71 



Table 9: Hadronic VP contributions to 10 °a M with FSR corrections included. Numbers within 
brackets refer to respectively statistical and systematic errors. Numbers within square brackets 
are the total uncertainties. 
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io 1( V 


Values 
Solution A 


(incl. r) 

Solution B 


Values ( 
Solution A 


excl. r) 

Solution B 


LO hadronic 


685.78 ±4.55 


687.72 ± 4.63 


682.82 ±4.66 


684.36 ±4.71 


HO hadronic 


-9.98 ± 0.04 exp ± 0.09 rad 


LBL 


10.5 ±2.6 


QED 


11 658 471.8096 ± 0.016 tot 


EW 


15.32 ±0.10 hadr ±0.15 Hi gg, 


Total Theor. 


11 659 173.43 ±5.25 


11 659 175.37 ±5.31 


11 659 170.47 ±5.34 


11 659 172.0 ±5.39 


Exper. Aver. 


11 659 208.9 ±6.3^ 


Aa M 


35.47 ±8.20 


33.53 ±8.24 


38.43 ± 8.26 


36.89 ±8.29 


Significance (na) 


4.33a 


4.07(7 


4.65(7 


4.45(7 



Table 10: The various contributions to 10 10 a M . Aa M = (a^) exp — (a^) th is given in units of 
10~ 10 and the last line displays its significance. 
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Figure 1: Best fits to e + e~ — > it + it~it cross sections for data sets in isolation. Left column 
displays fits of the CMD-2 data, right column displays fits of the SND data. Top shows the <f) 
region, bottom the u region. The plotted data are extracted from [|89l |901 (CMD-2) and [|92l 
(SND) for the region and from [|75] (CMD-2) and |ED (SND) for the u region. The empty 
circles (bottom right plot) are superimposed on the SND fit results and are not used in the fit 
displayed in this Figure. 
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Figure 2: Simultaneous fit of e + e~ — > n + n~7r° cross section on the region data from ||89ll90l 
(CMD-2) and flU (SND). 
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l-CH 


SND+CMD2 ^ (x 2 /N =176/1 13) 
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Figure 3: c\ — C2 values returned by fits. CMD2 u denotes the fit result of the data from 
lT75l SND u those from [91], ND+CMD the fit result to the merged data from [|9l and ll94l . 
CMD2 (ft indicates that only the merged data from [[88] |89] |90] have been used in the fit, SND </> 
corresponds to the fit of the data from [92 1 and SND+CMD2 <\> provides the (simultaneous) fit 
result of [88, 89[90l|92]. Finally, the last line shows the result for the selected data consisting 
of the sample reported in fl93l |9l |9]] [75] Ml Ml I2Q) . The vertical dotted line serves to show 
how the fits perform the averaging. 
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Figure 4: Simultaneous fit of the e + e — > n + n n data in the u and regions. Top figures 
show the case for the merged data from []9T] |92). Bottom figures display the fit results for 
CMD-2 data from [|75l[89ll90l. 
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Figure 5: Global fit of the e + e" — > 7r + 7r _ 7r° data . Top left enhances the u region, top right the 
4> region. The data superimposed are all fitted. Bottom plot shows the intermediate region; all 
plotted data are included in the fit procedure, except for the DM1 data set. The particular data 
sets used are described in the main text and in the captions to previous Figures. 
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Figure 6: Fit of the e + e — > KK data. Left side are K°K , right side K + K~ 
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Figure 7: Ratio of the e + e _ — > K°K and e + e~ — > K + K~~ cross sections normalized to the 
model ratio. Top panel displays the case for CMD-2 data, Bottom panel those for SND data. 
In the top panel, the residual experimental systematics band (2.3%) is figured by dashed lines. 
Correlated systematics between charged and neutral modes are expected to cancel out in the 
experimental ratios. 
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v's (GeV) 

Figure 8: Global Fit of the dipion spectrum in the decay of the r lepton. The data points are 
those from ALEPH [@0l, Belle @U and CLEO The inset magnifies the p peak region. 
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Figure 9: Global Fit of the pion form factor squared in e + e _ annihilations. The data points 
are those from CMD-2 [EH EH [77] and SND ED- 0ne has not plotted the so-called "old 
timelike" data also (mostly) collected at Novosibirsk. The inset magnifies the p peak region 
and the behavior at the p — to interference region. 
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Figure 10: Global Fit of the function H(s) = B W7T /NdN/ds in r decays. Top Figure shows the 
residuals as a function of s; downmost Figure shows the function (H fu(s) — Hdata(s)) / H /#(s) . 
The fitted region extends from threshold to 1.0 GeV/c, i.e. over the region where the behavior 
of the data sets from ALEPH [40], Belle [41Jand CLEO [|42ll reach some agreement. 
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Figure 11: Ratio of the transition amplitudes p° — 7 and — W ± , fpj/ f P w following from 
the global fit and neglecting loop corrections. This corresponds to the ratio shown in Tabled] 
and reproduced in Section [T21 Top Figure shows the real part as a function of s, bottom Figure 
the imaginary part. Uncertainties due to fit parameter errors are not given; the uncertainty band 
for f pl /fpw — 1 can estimated to a few percent. 
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Figure 12: Ratio of the couplings g W7r7r /g p7r7r as a function of ^fs, as coming from the global fit 
(this ratio is explicitly given in Section [13]). The vertical line locates the PDG mass of the u 
meson. The uncertainty band due to fit parameter errors is not shown. 
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Figure 13: A set of recent estimates of the muon anomalous magnetic moment together with 
the BNL average value JT1|2|. These are extracted from [21 (DHMZ10), El (JS11), ifTBll 
(HLMNT11) and |fT3l (DHea09). Our own results are figured by A and B for respectively 
solutions A and B. The statistical significance of the difference between the estimated and 
measured values of is displayed on the right side of the Figure for each of the reported 
analyses. 
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